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20.  ABSTRACT  (Continued) 

or  the  sphere  r singed  from  1100  to  85OO  ft/sec.  The  results  of  the  first 
year  of  investigation,  comprising  all  the  work  performed  at  China  Lake, 
have  been  documented  in  Report  DNA  3*+l6Z;  in  addition,  the  China  Lake 
performance  was  published  in  a  separate  report  NWC  TP  5700  (January  1975)* 


The  present  report  concentrates  on  the  activity  of  the  second  year 
of  operation.  Tests  were  performed  to  separately  measure  the  energy  con¬ 
tained  in  the  waves  generated  by  the  impact  in  the  sample,  and  special 
procedures  were  developed  and  executed  to  attempt  to  determine  the  surface 
area  associated  with  cracks  produced  as  the  result  of  the  collision. 

These  tests  were  performed  on  diorite  rather  than  on  shale,  since  the 
ejecta  formation  in  the  latter  was  found  to  be  negligible  for  the  range 
of  laboratory  impact  velocities  utilized.*  Photographic  measurements  of 
the  ejecta  and  the  determination  of  the  crater  topography  were  continued 
as  in  the  earlier  tests %  but  the  use  of  the  ballistic  pendulum  and  the 
scheme  for  the  evaluatiflPTif  the  spatial  distribution  of  the  ejecta  were 
abandoned^  Tests  on  shal^were  conducted  to  determine  the  energy  required 
for  given  depths  of  penetration . 

It  was  found  that,  over  a  ryige  of  six  decades  of  energy  of  the 
striker,  the  penetration  depth  in  shale  scaled  reasonably  well  with  both 
initial  momentum  and  initial  energy,  and  that  the  phenomenon  could  be 
satisfactorily  modeled  by  a  purely  viscous  target  material  that  neglects 
elastic  effects.  The  energy  of  comminution  for  the  harder  rock  diorite, 
consisting  of  the  difference  between  initial  striker  energy  and  the 
measured  terminal  energies  due  to  waves  in  the  target,  rebound  of  the 
projectile  and  ejecta  emission  were  found  to  be  of  the  order  of  JO  to 
85  percent.  In  view  of  the  extreme  uncertainties  as  to  what  cracks  in 
the  specimen  preexisted  and  what  fissures  were  directly  attributable  to 
the  impact,  it  was  only  possible  to  give  bounds  for  the  energy  required 
to  produce  unit  new  surface  in  the  phenomenon.  While  the  ejecta  surface 
area  was  determined  from  the  measured  size  distributions  and  a  geometric 
configuration  governed  by  a  single  linear  dimension,  a  complete  analytical 
representation  of  the  comminution  and  ejection  process  will  require  con¬ 
siderable  further  study  and  experimentation. 


si 


1 

j 


_ UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  THIS  PAGEfWltn  0»«  Entind) 


TABI£  OF  CONTENTS 


Page 


Section  I  Introduction  5 

Section  II  Summary  of  First  Year  Activity  7 

Section  III  Liaison  with  Stanford  Research  Institute  ^ 

Section  IV  Work  performed  during  the  Second  Contract  Year  ^ 

IV. A  Sample  Preparation  ^5 

IV.  B  Experimental  Apparatus  and  Procedure  23 

IV. C  Results  and  Discussion  .  33 

1.  Diorite  33 

a.  Velocity  Measurements  33 

b.  Ejecta  36 

c.  Wavj  Energy  53 

d.  Crater  56 

e.  Cracks  56 

f.  Relation  of  Results  to  Comminution  Theories  69 

2 .  Shale  72 

a.  Ejecta  72 

b.  Dynamic  Penetration  74 

c.  Quasi-static  Tests  77 

d.  Correlation  of  Penetration  Data  8l 

Section  V  Summary 

List  of  References  ®7 

Appendix  A  Computer  Program  for  Crack  Surface  Area  Evaluation  39 


1 


Figure 


LIST  OF  ILLUSTRATIONS 


Title 


Composite  Specimen  with  Embedded  Gages 


Assembly  of  Composite  Diorite  Sample 


Schematic  of  Experimental  Arrangement 


Photograph  of  Experimental  Arrangement 


Typical  Record  of  an  Initial  Velocity  Measurement 


Strain  Gage  Histories  for  Run  D180 


Montage  of  Typical  Cross-section  of  Diorite  Sample 
after  Impact 


Sectioning  Scheme  for  Diorite  Target  Plates 


Typical  High-speed  Photograph  of  Diorite  Ejecta 


Assumed  Ejecta  Configuration  and  Velocity  Distribution 


Typical  Ejecta  as  Viewed  by  Quantimet  720  Unit 


Histogram  of  Particle  Fraction  versus  Size  for  Run  D180 
Assumed  Spherical  Particle  Shape 


Probabilistic  Ejecta  Distribution  for  Run  D180 


Energy  Input  versus  Mean  Particle  Size  for  Three  Particle 
Shapes 


Comparison  of  Measured  and  Predicted  Tangential  Strain 
Histories  for  Run  D180 


Typical  Contour  Map  of  a  Diorite  Crater 


Typical  CALCOMP  Plot  of  the  Cross-sectional  Crack 
Distribution 


Typical  Photographs  of  Dyed  Samples 


High-speed  Photograph  of  a  Typical  Shale  Test 


Section  of  Shale  Sample  Struck  by  a  1/4  inch  Diameter  Steel 
Sphere  at  a  Velocity  of  5300  ft/sec.  Crater  filled  with  a 
catalyzed  methylmethacrylate  glue  prior  to  sectioning. 


Til 


2 


Figure 

Title 

Page 

20 

Penetration  Depth  as  a  Function  of  Initial  Projectile 
Momentum  for  Shale 

82 

21 

Square  of  Penetration  Depth  as  a  Function  of  Initial 
Projectile  Energy  for  Shale 

83 

3 


2to*!Ssifc8J ass 


rjissm"±:~ 


'i-mm Mmmamsm*  ?  ;  z  • 


LIST  OF  TABLES 

Table 

Title 

Page 

1 

Mechanical  Properties  of  Shale 

16 

2 

Mechanical  Properties  of  Diorite 

17 

3 

Impact  Plate  Thickness  and  Strain  Gage  Locations 

21 

4 

Semi-Conductor  Strain  Gage  Specifications 

22 

5 

Ejecta  Cloud  Characteristics 

37 

6a 

Ejecta  Calculations  for  Run  D180. 

Disk  with  height  =0.3  radius 

43 

6b 

Ejecta  Calculations  for  Run  D180. 

Disk  with  height  =0.5  radius 

44 

6c 

Ejecta  Calculations  for  Run  D180.  Sphere 

45 

7 

Mean  and  Standard  Deviation  of  Ejecta  Surface 

Areas  for  Three  Particle  Shapes 

48 

8 

Diorite  Test  Results 

50 

9 

Summary  of  Wave  Energy  Calculations 

57 

10 

Crack  Surface  Areas  Based  on  Various  Assumptions 

67 

11 

Comminution  Energy  per  Unit  of  New  Surface  Area 

70 

12 

Comminution  Constants  Based  on  Quasi-static  Theories 

71 

13 

Dynamic  Shale  Penetration  Data  for  a  ^-in.  Diameter 

Steel  Sphere 

75 

14 

Quasi-static  Penetration  Tests  on  Shale 

78 

4 


SECTION  I 
INTRODUCTION 

The  state  of  the  art  in  the  field  of  fragmentation  or  comminution  of 
rock  is  not  highly  developed,  yet  this  process  is  of  extreme  importance  in 
a  variety  of  military  and  civilian  applications.  Furthermore,  the  phenomena 
incident  to  comminution  produced  by  impact  are  even  less  well  understood 
than  those  associated  with  quasi-static,  more  or  less  steady  state  events 
such  as  occur  in  mechanical  size  reduction  procedures,  yet  the  former  consti¬ 
tute  an  area  of  considerable  interest  to  those  concerned  with  the  effects  of 
weapons.  In  particular,  the  impact  of  a  projectile  or  the  transient  loadings 
caused  by  an  offensive  device  on  the  surface  of  a  consolidated  geological 
material  produce  extremely  complicated  phenomena.  The  present  document 
reports  an  investigation  that  was  designed  to  provide  some  information  on 
the  mechanical  aspects  of  this  problem  believed  to  be  of  importance  in  an 
understanding  of  the  dynamic  comminution  process. 

Comminution  as  a  field  of  study  has  a  history  of  more  than  one  century, 

but  has  been  almost  exclusively  directed  towards  industrial  applications. 

Two  early  theories  pertaining  to  the  energy  required  for  the  sisie  reduction 

process  were  proposed  by  Rittinger^  and  Kick^;  a  third  hypothesis  was 

f3] 

advanced  70  years  later  by  BondL  ,  and  a  more  recent  version  of  a  more 

[4] 

inclusive  nature  was  presented  by  Walker  .  These  expressions  are  based 
principally  on  dimensional  analyses  and  do  not  pretend  to  portray  an  appro¬ 
priate  phenomenological  model  of  the  process.  The  vast  reservoir  of  experi¬ 
mental  information  derived  in  the  interim  has  been  concerned  with  energy 
requirements  for  crushing  and  grinding  operations  and  the  resultant  particle 
size  distributions  for  a  variety  of  brittle  materials Virtually 
none  of  these  studies  relate  to  the  problem  of  rock  disintegration  under 
conditions  of  impact,  and  it  cannot  be  assumed  a  priority  that  such  quasi¬ 
static  results  can  be  applied  to  projectile  impact  on  geological  materials. 

Such  an  event  involves  a  complicated  sequence  of  mechanical  phenomena 
that  predominate  over  thermal  effects,  in  contrast  to  the  case  of  nuclear 
explosions  where  thermal  and  radiation  energies  are  highly  significant. 
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The  impact  case  consists  of  a  local  dynamic  contact  problem  which  results  in 
wave  propagation  in  the  target  that  produces  motion,  strain,  cracking,  ejec 
tion  of  particles,  compaction  and  crater  formation.  It  may  also  produce 
plastic  deformation  and/or  shattering  of  the  striker  when  certain  combinations 
of  high-velocity  impact  and  material  hardness  prevail;  the  former  did  not 
occur  to  any  measurable  degree  and  the  latter  situation  was  found  only  in  one 
special  sequence  of  the  present  series  of  experiments.  The  investigation 
described  addresses  itself  to  the  measurement,  data  analysis  and  explanation 
of  these  various  aspects. 

This  document  constitutes  the  final  technical  report  for  the  work  per¬ 
formed  under  contract  No.  DNA001-73-C-0226  between  the  University  of  Cali¬ 
fornia,  Berkeley,  and  the  Defense  Nuclear  Agency  who  supported  the  activity 
under  NWE  D,  Subtask  Code  Y990AXSA001,  Work  Unit  Code  42.  It  completes  the 
work  described  in  document  DNA  3416 Z,  an  interim  report  for  the  period  15 
May  1973  -  14  May  1974,  and  covers  the  interval  up  to  30  November  1975^. 

The  essential  features  of  the  research  involved  the  following  topics: 

1.  The  examination  of  the  process  of  comminution  in  rocks  by  experi¬ 
mentally  determining  the  mechanisms  relating  to  the  fragmentation  process 
when  a  striker  impinges  the  plane  surface  of  a  block  of  this  material.  In 
part,  this  was  to  be  accomplished  through  the  use  of  high-speed  photography 
that  would  provide  a  measure  of  the  ejecta  velocity  and  angle  of  ejection  as 
well  as  the  projectile  velocities.  The  particle  size  distribution  was  deter¬ 
mined  by  analysis  of  the  collected  ejecta  and  its  spatial  distribution  was 
also  ascertained. 

2.  An  energy  balance  of  the  overall  phenomenon  was  performed  in  order 
to  attempt  to  obtain  a  value  of  the  energy  required  to  produce  unit  new  sur¬ 
face.  To  execute  this  calculation,  it  was  necessary  to  measure  the  kinetic 
and  strain  energy  transmitted  to  the  target  by  wave  processes  and  to  assess 
the  magnitude  of  the  new  surfaces  generated.  This  involved  the  employment  of 
a  simple  model  of  the  ejecta  geometry  as  characterized  by  a  single  spatial 
parameter  in  conjunction  with  the  measured  particle  size  distribution.  It 
was  also  possible  to  identify  those  aspects  of  the  phenomenon  exhibiting 
very  small  energy  components  that  could  be  neglected  in  the  calculations. 
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3.  Within  the  framework  of  the  overall  effort  on  this  problem,  the 
possibility  of  scaling  the  event  in  terms  of  frequently-employed  mechanical 
variables  was  examined  for  wide  variations  in  velocity  and  geometries  of 
the  strikers. 

While  the  objective  of  the  task  was  initially  directed  towards  the 
examination  of  the  behavior  pattern  in  soft  rocks,  considerable  effort  was 
also  expended  on  an  analogous  study  of  hard  rock  response  to  impact.  Tests 
involving  the  first  type  of  target  required  the  measurement  of  projectile 
penetration  as  well  as  ejecta,  although  the  latter  was  found  to  be  very 
small  in  a  number  of  experiments  implying  that  the  major  mechanism  of  pene¬ 
tration  was  due  to  material  compaction;  some  quasi-static  penetration  studies 
were  also  performed.  In  constrast,  the  ejecta  formed  in  the  hard  rocks 
produced  sizable  shallow  craters  whose  topography,  including  volume,  were 
ascertained. 

The  work  on  this  project  involved  a  collaborative  effort  with  the  Naval 

Weapons  Center,  China  Lake,  California  that  was  completely  reported  in  the 

interim  document  and  will  not  be  repeated  here.  Close  cooperation  was  also 

maintained  with  the  Geophysics  Group  at  the  Stanford  Research  Institute, 

Menlo  Park,  California,  as  requested  in  the  contract  for  the  second  year. 

A  series  of  experiments  involving  novaculite  was  performed  for  this  group 

in  the  laboratories  at  Berkeley  whose  results  are  presented  in  the  final 

r  9] 

report'  .  This  activity  required  considerable  modification  of  the  equipment 
and  procedure  that  will  be  briefly  described  in  the  sequel. 

SECTION  II 

SUMMARY  OF  FIRST  YEAR  ACTIVITY 

The  investigation  during  the  first  year  of  the  contract  was  performed 
both  in  the  laboratories  at  the  University  of  California,  Berkeley,  and  at 
the  Naval  Weapons  Center,  China  Lake,  California.  While  a  full  and  complete 
account  of  this  activity  is  described  in  Reference  [8],  a  summary  of  this 
work  featuring  the  highlights  and  essential  results  will  be  presented  in 
the  sequel. 

The  principal  experimental  apparatus  for  the  comminution  studies  in 
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the  Berkeley  laboratories  consisted  of  a  ballistic  pendulum  with  attached 
target  that  was  struck  by  a  1  /?.  inch  diameter  steel  sphere  projectile  fired 
from  a  pneumatic  propulsion  device  at  velocities  ranging  from  about  450  to 
700  ft/sec.  The  velocity  of  the  projectile  was  measured  by  the  signal  pro¬ 
duced  from  the  interruption  of  two  light  beams,  6  inches  apart,  traversing 
two  slots  in  the  gun  barrel  near  the  muzzle  end.  Passage  of  the  striker  upon 
exit  from  the  barrel  through  the  beam  of  a  laser  also  provided  a  signal  which, 
upon  proper  delay,  activated  the  discharge  of  a  condenser  bank  into  a  flash 
unit  serving  as  the  illumination  for  the  photographic  recording  of  the  event 
by  a  high-speed  Beckman-Whitley  framing  camera.  The  targets  consisted  of 
shale,  limestone  or  diorite  rock  samples  with  a  diameter  of  5-1/2  inches  and 
a  thickness  of  1-1/2  inches,  the  first  being  encased  in  a  matrix  of  plaster 
of  Paris,  while  the  latter  was  unconfined;  both  types  of  samples  had  a  short 
aluminum  cylinder  glued  to  the  center  of  their  top  surfaces  for  attachment  to 
the  pendulum.  The  swing  of  the  latter  was  measured  by  a  potentiometer  and 
recorded  on  an  oscilloscope. 

The  ejecta  produced  by  the  impact  formed  a  distinctive  cloud,  a  portion 
of  which  was  trapped  on  a  10-1/8  inch  diameter  lucite  plate  coated  with  a 
polymeric  agar  to  a  depth  of  1/4  inch,  located  approximately  2  inches  from 
the  initial  pendulum  position.  Particles  not  caught  on  this  plate  were 
arrested  by  a  plastic  tube  with  an  inside  diameter  of  11-3/4  inches  and  a 
length  of  21-1/4  inches  that  completely  enveloped  the  rock  sample  and  plate. 

The  photographic  data,  obtained  on  79  frames  each  with  dimensions  of 
about  1  cm  x  2.5  cm  at  rates  ranging  from  50,000  to  190,000  frames  per  second 
permitted  the  determination  of  the  initial  and  final  striker  velocity,  the 
contact  history  of  the  ball,  and  the  characteristics  of  the  debris  generated. 
The  ejecta  formed  a  cloud  of  fine  particles  as  well  as  some  larger  fragments 
that  could  be  individually  identified  in  a  sequence  of  frames.  The  speed  of 
the  ejecta  was  calculated  in  the  plane  normal  to  the  camera,  that  for  the 
front  of  the  rock  cloud  being  taken  as  the  average  of  the  values  along  three 
rays  with  due  consideration  for  the  shape  of  this  cloud.  The  softer  rocks, 
shale  and  limestone,  both  exhibited  a  cylindrical  ejecta  pattern,  while  that 
of  the  diorite  formed  a  conical  shell.  The  velocity  of  the  larger  fragments 
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was  also  determined  on  the  basis  of  uniform  rectilinear  motion,  although  this 
is  only  a  crude  approximation  in  view  of  the  known  particle  deceleration  and 
some  observed  rotation  of  the  larger  fragments. 

Results  from  tne  film  also  permitted  the  calculation  of  the  initial  and 
rebound  energy  of  the  striker  together  with  an  estimate  of  the  kinetic  energy 
of  the  ejecta.  The  particles  entrapped  on  the  catcher  plate  were  analyzed 
with  respect  to  spatial  distribution  by  means  of  an  electronic  scanner,  which 
sampled  specified  plate  regions  with  respect  to  particle  size,  area,  peri¬ 
meter,  and  a  shape  factor  relative  to  a  circular  configuration.  Particles 
not  so  caught,  but  collected  in  the  plastic  tube  were  weighed  and  also 
subjected  to  analysis  by  the  scanning  unit. 

An  energy  balance  was  executed  on  the  basis  of  the  relation 

bi  bf  r  cl  pend  com  fr  pi  o  v  1 

where  the  first  term  is  the  initial  striker  energy,  and  the  terms  on  the 
right-hand  side  represent,  respectively:  the  terminal  striker  energy;  that 
of  the  larger,  separately  identified  rock  fragments;  the  kinetic  energy  of 
the  cloud;  that  of  the  pendulum;  the  energy  of  comminution  (i.e.,  separation 
energy  of  the  ejecta),  the  energy  of  crack  formation  in  the  target;  that  of 
plastic  target  and  projectile  deformation;  and  other  energies,  principally 
that  of  elastic  waves  in  hard  rock  targets.  In  the  initial  calculations, 


the  energy  of  plastic  deformation  and  wave  energy  were  neglected.  In  none 
of  the  first  year  tests  was  any  evidence  found  to  indicate  the  presence  of 
projectile  deformation.  However,  estimates  based  on  strain  measurements  on 
the  target  surface  showed  that  the  wave  energy  in  hard  rocks  was  not 


negligible.  Furthermore,  although  no  plastic  deformation  of  the  hard  rock 
samples  was  noted,  a  great  deal  of  compaction  was  found  to  exist  in  the  shale 
targets.  With  the  omission  of  the  above  energies,  it  was  stated  that  99 
percent  of  the  initial  kinetic  energy  of  the  striker  was  converted  into 
comminution  and  fracture  for  the  shale,  although  this  figure  in  reality  also 
includes  the  very  sizeable  energy  of  deformation  component,  dn  the  other 
hand,  for  limestone  and  diorite  where  only  negligible  compaction  was  observed, 
the  comminution  and  fracture  energy  was  found  to  be  94  percent  and  80  percent, 
respectively,  based  on  the  same  method  of  calculation. 
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The  significant  amount  of  material  compaction  in  the  shale  samples  was 
evidenced  by  the  fact  that  circular  lines  drawn  initially  around  the  impact 
point  were  subsequently  found  at  the  bottom  of  the  nearly  uniform  depression. 
The  crater  topography  was  measured  by  stereotactic  means  over  a  discrete  net 
of  points  using  a  feeler  gage.  The  data  obtained  in  this  manner  constituted 
the  input  to  a  computer  program  that  calculates  both  hole  volume  and  a  con¬ 
tour  map.  The  crater  volume  was  also  determined  by  means  of  fluid  o'*  solid 
filler  materials  for  compact  and  porous  rocks,  respectively,  and  in  “tion, 
by  an  estimate  based  on  micrometric  measurements  of  the  maximum  and  mi'.,  '  l.  i 
dimensions  of  the  crater  near  the  surface  and  the  maximum  depth.  The  shape 
of  the  crater  in  shale  was  nearly  hemispherical,  with  a  circular  surface 
contour,  while  that  for  diorite  was  conically  shaped;  the  volume  for  the 
former  material  was  about  one  order  of  magnitude  larger  than  for  the  harder 
rock  in  the  case  of  identical  impact  situations. 

The  surface  area  distribution  for  fields  examined  at  radii  of  5,  50,  100 
and  105  mm  from  the  center  of  the  catcher  plate  could  generally  be  represented 
by  a  triangular  (or,  alternatively,  by  a  Gaussian)  distribution  on  a  loga¬ 
rithmic  scale  for  angles  of  ejection  of  0,  45,  90,  135,  180,  225,  270  and 

2 

315  degrees.  For  a  range  of  surface  areas  from  about  20  to  60,000  y  the  naxi- 
nun  nutter  of  particles  involved  was  140,  80,  80,  50,  40,  40,  70,  and  210, 
respectively  (See  Figs.  12a-12h,  Reference  [8]).  Most  of  the  particles  were 
trapped  in  the  region  of  the  largest  plate  radius  rather  than  in  the  central 
region;  thus,  many  shale  fragments  were  not  caught  by  the  arresting  plate 
as  was  also  indicated  by  the  substantial  amount  of  debris  remaining  in  the 
tube. 

The  zonal  surface  area  of  the  particles  for  the  shale  was  found  to  vary 
linearly  with  radius  of  the  catcher  plate  except  when  the  latter  was  moved 
very  close  to  the  target;  in  this  case,  the  peak  value  of  the  annular  surface 
area  did  not  occur  at,  albeit  near,  the  plate  edge.  However,  based  on  an  extra¬ 
polation  of  the  data  from  the  ratio  of  the  area  scanned  by  the  Quantimet 
unit  to  that  of  the  entire  catcher  plate,  about  18  percent  of  the  crater  mass 
of  the  shale  sample  was  arrested  by  the  plate,  with  only  an  additional  10 
percent  of  the  mass  recovered  as  fragments  in  the  enveloping  plastic  tube. 
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While  a  large  number  of  dust  particles  undoubtedly  escaped  from  the  system, 
the  major  reason  for  this  extreme  discrepancy  is  undoubtedly  due  to  the 
compaction  of  the  target. 

The  work  performed  at  China  Lake,  also  detailed  in  Reference  [8],  con¬ 
sisted  of  three  separate  types  of  tests,  including  small-scale  laboratory 
investigations  involving  1/4  inch  diameter  spherical  steel  projectiles, 
firings  of  20  and  40  mm  diameter  projectiles  from  a  cannon  against  the  verti¬ 
cal  shale  cliff  of  the  open  pit  mine  of  the  American  Borax  Company  at  Boron, 
California,  and  a  single  firing  of  a  6  inch  diameter  projectile  into  a  reddish 
silt  whose  geological  properties  were  believed  to  closely  resemble  those  of 
shale. 

The  small-scale  firings  were  conducted  principally  into  shale  disk  tar¬ 
gets,  either  plaster-of-Paris  encased  or  enveloped  by  tape,  with  diameters  of 
5.5  inches  and  thicknesses  ranging  from  1-1/4  inch  to  5.7  inches  at  velocities 
up  to  8500  ft/sec.  Most  of  the  specimens  were  backed  by  a  metallic  disk  of 
identical  diameter  that  prevented  rearward  spalling .  The  samples  were  mounted 
either  on  a  ballistic  pendulum  or  else  on  a  rigid  test  stand,  the  excursion 
of  the  former  being  measured  by  a  precision  potentiometer  attached  to  the  axis 
of  rotation.  The  phenomenon  was  photographed  by  means  of  a  six-frame  Beckman- 
Whitley  Kerr  cell  type  of  camera,  and  the  debris  liberated  by  the  impact  was 
collected  in  a  polyethylene  bag  draped  over  the  target  configuration. 

The  impact  on  the  shale  produced  a  cylindrical  hole  carved  out  by  the 
projectile  almost  exclusively  as  the  result  of  material  compaction  as  well 
as  a  cone-shr ped  or  dish-like  depression  near  the  surface  apparently  formed 
by  ejection  of  particles.  By  superposing  a  metal  plate  with  a  central  hole 
to  permit  passage  of  the  projectile  on  the  sample,  the  surface  dishing  was 
virtually  eliminated.  The  ejecta  collected  under  these  circumstances  amounted 
to  20-25  percent  of  the  total  hole  volume,  independent  of  initial  velocity; 
the  remainder  of  the  volume  had  to  be  attributed  to  compaction.  On  the  other 
hand,  only  dishing  occurred  in  the  harder  and  more  brittle  rocks  such  as 
diorite  and  spessartite,  with  some  propagation  of  cracks  produced  by  the 
impact;  compaction,  plastic  deformation  or  distortion  of  both  target  and 
striker  are  so  small  as  to  be  neglegible. 
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Sectioning  of  the  damaged  targets  was  performed  after  filling  the  hole 
with  a  clear  epoxy  mixture  to  avoid  further  breakage.  This  filler  also 
clearly  delineated  the  region  of  compaction,  which  consisted  of  a  bell-shaped 
zone  around  the  hole  with  its  tip  approximately  one  radius  ahead  of  the  ter¬ 
minal  projectile  position. 

A  total  of  46  rounds  was  fired  in  this  sequence  of  experiments,  in¬ 
cluding  preliminary  tests  involving  diorite  and  spessartite  targets.  The 
crater  depth  was  found  to  vary  linearly  and  the  crater  volume  as  the  square 
of  the  initial  velocity,  with  large  deviations  occurring  only  upon  shattering 
of  the  projectile;  the  latter  event  occurred  only  at  extreme  velocities.  The 
volume  of  the  surface  dish  was  found  to  vary  between  30  and  80  percent  of  the 
total  crater  volume  without  any  apparent  relation  to  initial  velocity;  in  all 
cases,  an  extensive  network  of  cracks  was  produced  that  frequently  extended 
to  the  lateral  specimen  surfaces. 

The  fragment  velocity  was  found  to  increase  with  rising  initial  projec¬ 
tile  speed,  and  the  angle  of  ejection  was  noted  to  decrease  correspondingly 
from  a  magnitude  of  about  60°  in  the  range  from  2000  to  3000  ft/sec  to  an 
angle  of  35°  at  the  maximum  velocity  of  about  8200  ft/sec.  In  addition  to 
the  symmetric  particle  cloud  forming  an  outer  conical  shell  that  was  subse¬ 
quently  transformed  into  a  double  envelope,  a  central  initially  circular 
region  of  particles  was  photographed  that  was  believed  to  be  the  result  of 
high  turbulence  inducing  vortex  formation.  The  sieving  of  the  debris  by  a 
series  of  meshes  did  not  result  in  a  Gaussian  distribution. 

A  total  of  18  rounds  was  fired  from  tank  guns  using  20  mm  rifled  and  40 
mm  rifled  and  unrifled  projectiles  with  standard  ogival  noses  against  fresh 
vertical  shale  surfaces  of  the  Boron  open  pit  mine  at  initial  velocities 
ranging  from  2000  to  2900  ft/sec  as  determined  by  the  piercing  of  electrically 
wired  screens.  This  provided  an  initial  kinetic  energy  ranging  from  0.30- 
2.33  x  10  5  ft-lbs  compared  to  a  value  of  0.000012  x  10  5  ft-lbs  for  the  1/4 
inch  diameter  steel  spheres  fired  at  185  ft/sec.  The  ph  :.nenon  was  photo¬ 
graphed  by  two  Fastax  caneras  operating  at  nominal  framing  rates  of  3000  per 
second,  indicating  a  conical  cloud  angle  of  about  60°  and  maximal  fragment 
velocities  of  about  700  ft/sec.  After  clearing  the  base  of  the  cliff,  the 
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debris  ejected  by  the  firing  of  the  round  was  collected  in  large  plastic 
sheets  and  weighed;  however,  a  portion  of  this  material  from  remote  regions 
of  the  wall  was  believed  to  have  been  knocked  loose  by  either  the  sound  or 
shock  waves.  While,  for  the  40  mm  projectiles  at  full  loads,  considerable 
scatter  was  found  in  the  penetration  depth  and  the  type  of  trajectory  (a 
number  of  projectiles  veered  from  their  path  or  began  to  tumble),  consistent 
results  were  obtained  for  40  mm  projectiles  with  reduced  charges  and  for  the 
20  mm  projectile.  Surface  dish  volume  data  for  all  shots  were  also  incon¬ 
sistent;  this  might  be  expected  in  view  of  the  variable  and  uncontrollable 
surface  conditions  of  the  various  parts  of  the  cliff  where  shots  were 
attempted.  Significant  plastic  deformation  of  the  cliff  was  found  for  all 
rounds  at  the  entry  into  the  bore  hole  without,  however,  leading  to  changes 
in  density  of  the  substance. 

A  fair  correlation  exists  between  the  initial  kinetic  energy  of  these 
larger  shots  and  the  weight  of  shale  from  the  crater;  here,  most  of  the 
material  so  removed  came  from  the  dish,  while  a  larger  proportion  was  dis¬ 
placed  from  the  hole  in  the  small-caliber  tests.  This  is  partly  due  to 
lateral  confinement  of  the  model  and  the  absence  of  loose  or  projected 
sections  in  the  machined  samples. 

Sieving  of  the  recovered  debris  also  did  not  yield  a  Gaussian  size  or 
weight  distribution;  undoubtedly,  secondary  fragmentation  of  the  material 
occurred  to  some  extent  upon  contact  with  the  ground.  Only  a  small  fraction 
of  the  weight  of  the  material  removed  from  the  crater  based  on  calculations 
involving  the  hole  volume  and  the  initial  density  was  recovered;  the  remainder 
was  ejected  over  an  area  much  wider  than  the  ground  cloth,  and,  most  impor¬ 
tantly,  sizeable  compression  and  plastic  deformation  of  the  cliff  ensued. 

A  single  firing  of  a  6  inch  diameter  ogival-nosed  steel  shell  weighing 

3 

95.6  lb.  into  a  cohesive  reddish  silt  having  a  density  of  1.52  g/cm  was 
executed  at  a  velocity  of  1340  ft/sec.  The  trajectory  was  not  linear,  with 
a  total  travel  of  about  11  feet.  The  dished  region  formed  a  cone  with  a 
surface  diameter  of  3  feet  and  a  depth  of  6  inches;  the  initial  diameter  of 
the  penetrating  hole  was  identical  to  that  of  the  projectile.  The  total 

3 

crater  volume  is  thus  5615  in  with  a  weight  of  308  lb  for  an  initial 
kinetic  energy  of  2.68  x  10  ^  ft-lb. 
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SECTION  III 

LIAISON  WITH  STANFORD  RESEARCH  INSTITUTE 

As  requested  in  the  contract,  close  liaison  was  maintained  throughout 
the  year  with  the  Geophysics  Group  of  Stanford  Research  Institute,  particularly 
Drs.  D.  Curran  and  D.  Shockey.  At  their  request,  a  series  of  tests  were  per¬ 
formed  on  the  hard  rock  novaculite  using  the  pneumatic  gun  and  high-speed 
camera  of  the  Berkeley  laboratory.  While  the  samples  were  furnished  by  the 
Stanford  group,  the  experiments  required  considerable  modification  of  the 
test  arrangement,  particularly  with  respect  to  the  ejecta  recovery  unit. 
Additional  changes  were  required  to  accomodate  a  device  placed  on  the  muzzle 
end  of  the  gun  that  insured  normal  impact  of  the  projectile  with  the  specimen 
and  simultaneously  determined  the  exit  velocity.  The  tool  steel  projectiles 
for  these  experiments  shattered  upon  impact  due  to  the  high  velocity  and 
their  brittle  nature.  Two  of  these  tests  were  performed  and  the  data  was 
transmitted  to  the  requesters.  The  results  have  been  described  in  the  final 
report  of  the  group.  Reference  [9]. 

SECTION  IV 

WORK  PERFORMED  DURING  THE  SECOND  CONTRACT  YEAR 

The  principal  activity  during  the  second  year  of  operation  of  the  con¬ 
tract  was  the  determination  of  the  energy  required  to  gc  ierate  unit  new 
surface  in  the  case  of  impact  on  hard  rocks,  the  evaluation  of  the  penetra¬ 
tion  of  projectiles  in  the  case  of  soft  rocks,  and  an  attempt  to  compare  the 
data  obtained  during  the  entire  course  of  the  investigation  with  both  analy¬ 
tical  expressions  developed  for  fragmentation  (primarily  intended  to  describe 
grinding  processes)  and  with  parameters  scaling  the  impact  intensity.  This 
involved  several  relatively  independent  physical  aspects  of  the  problem: 

(a)  The  determination  of  the  energy  contained  in  the  elastic  wave  propagating 
through  the  specimen,  (b)  The  evaluation  of  the  surface  area  associated  with 
the  ejecta,  continuing  the  endeavors  of  the  first  year  of  operation,  (c)  The 
estimation  of  the  new  surface  generated  by  the  cracks  produced  in  the  specimen 
due  to  the  impact,  (d)  The  crater  morphology,  and  for  the  soft  rock,  the 
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penetration  depth.  Examination  of  the  projectiles  subsequent  to  contact 
indicated  no  permanent  deformation  for  any  of  these  tests.  In  the  sequel, 
the  details  of  this  effort  will  be  described. 

IV. A.  SAMPLE  PREPARATION 

The  rock  specimens  employed  consisted  of  samples  of  green  shale  and 
diorite,  5-1/2  inch  in  diameter,  both  cored  by  my  means  of  a  standard  rock 
drill.  The  shale  was  obtained  from  the  site  of  the  open  pit  mine  of  the 
U. S.  Borax  and  Chemical  Corporation,  Boron,  California;  the  geology  and 
details  of  sample  preparation  are  give  in  Reference  [8].  The  samples  were 
enclosed  in  plaster  of  Paris  molds  with  a  6.5  inch  square  cross  section;  the 
specimens  were  cut  with  a  water-cooled  diamond  saw  to  a  thickness  of  1  inch. 
The  perimeter  of  the  samples  was  wrapped  with  friction  tape  to  prevent  lateral 
cracking  of  the  specimens  and  mold  during  testing.  Table  1  lists  pertinent 
mechanical  properties  that  have  been  derived  for  this  material. 

The  diorite  acquisition  site,  mineralogy  and  sample  preparation  were 
also  described  in  Reference  [8];  the  same  procedure,  but  without  encasement, 
was  employed  to  obtain  1  inch  thick  samples  with  the  coring  and  cutting  at 
feeding  rates  of  a  few  inches  per  hour.  Table  2  presents  the  mechanical 
properties  of  this  material;  additional  information  on  the  mineralogy  and 
properties  of  this  substance  are  contained  in  References  [10]  and  [11]. 

The  determination  of  the  energy  contained  in  elastic  waves  transmitted 
to  the  sample  required  the  construction  of  a  reusable  specimen  with  strain 
gages  embedded  in  the  interior.  This  could  not  be  accomplished  for  the 
shale  due  to  its  fragility.  As  shown  in  Figure  1,  this  composite  consisted 
of  a  very  thin  frontal  plate  of  diorite  and  a  somewhat  thicker  back-up  plate 
of  this  material,  both  of  the  same  diameter,  5.5  inch,  as  the  gaged  1  inch 
thick  specimen.  The  latter  was  cut  on  a  diametral  plane  with  a  0.08  inch 
thick  water-cooled  diamond  saw,  lapped  to  4/0  grain  size  using  standard 
procedures,  and  three  radial  strain  gages  were  then  cemented  along  the  axis 
of  symmetry  as  indicated.  No.  38  enamel-coated  leads  with  a  diameter  of 
0.004  inch  and  a  resistance  of  0.660  ohm/ foot  were  attached  to  the  gage  wires 
and  brought  out  along  a  diameter  of  the  section.  The  two  halves  of  the 
specimen  were  then  glued  together  with  3M  structural  epoxy,  forming  a  joint 
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0.1  inch  thick  that  was  cured  for  at  least  one  day  and  whose  rough  edges  were 
smoothed  with  increasing  grades  of  fineness  of  emery  cloth.  One  tangential 
gage  was  then  glued  to  each  of  the  top  and  bottom  surfaces  at  the  center, 
as  shown  in  Figure  1,  and  lead  wires  were  again  attached  and  brought  to  the 
outside  of  the  composite  in  a  similar  manner. 

Several  other  1  inch  x  5.5  inch  diameter  diorite  disks  were  further  cut 
to  various  thicknesses  with  a  1/8  inch  thick  oil-cooled  diamond  saw  at  a 
rate  of  a  few  inches  per  hour.  Two  of  these  were  the  0.035  inch  thick  target 
face  plate  and  the  1/4  inch  thick  back-up  plate  which  were  permanently  bonded 
to  the  front  and  back  of  the  gaged  section,  respectively,  after  both  sides  had 
been  lapped  as  described  above.  A  nominal  3/4  inch  thick  diorite  transmission 
plate  was  then  wrung  on-  with  grease  to  the  rear  of  the  back-up  plate,  and  a 
number  of  impact  plates  ranging  in  thickness  from  0.205  inch  to  0.434  inch, 
that  essentially  constituted  the  sample  undergoing  comminution,  was  also  wrung 
on  with  grease  to  the  front  of  the  target  face  plate.  Finally,  this  diorite 
composite  was  wrung  with  grease  onto  an  11  inch  diameter  aluminum  back-up 
plate  consisting  of  2  disks  with  a  total  thickness  of  4.5  inches  to  serve  as 
a  wave  trap.  To  insure  intimate  contact  of  the  greased  faces,  the  unit  sat 
under  its  own  weight  for  some  time  prior  to  testing,  so  as  to  permit  the 
escape  of  air  bubbles  and  excess  grease  from  the  joints.  The  complete  assem¬ 
bly  is  shown  in  Figure  2;  Table  3  gives  a  description  of  the  specimen  proper¬ 
ties  and  the  distance  of  the  five  strain  gages  from  the  impact  point  which 
was  located  at  the  center  of  the  front  surface  of  the  impact  plate. 

The  five  strain  gages,  denoted  as  1R,  2R  and  3R  for  the  radial  trans¬ 
ducers  and  IT  and  2T  for  the  tangential  ones,  consisted  of  SPB-3-12-12  semi¬ 
conductor  gages  manufactured  by  Baldwin-Lima-Hamilton  with  a  nominal  gage 
length  of  1/8  inch,  an  unstrained  resistance  of  128.2  ohms  +  1  percent  and  a 
gage  factor  of  117.3,  according  to  the  manufacturer.  The  gages  were  cemented 
with  Baldwin-Lima-Hamilton  gage  epoxy.  Table  4  presents  the  measured  circuit 
resistance  for  each  transducer  at  two  times  during  the  period  of  investigation. 
Each  gage  was  incorporated  in  a  potentiometric  circuit  activated  by  a  24-volt 
battery  with  a  2120  ohm  ballast  resistor  that  reduced  the  voltage  across  each 
gage  and  thus  prevented  their  overheating.  Each  gage  was  calibrated  after 
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Assembly  of  Composite  Diroite  Sample 
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TABLE  4.  SEMI-CONDUCTOR  STRAIN  GAGE  SPECIFICATIONS 
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supplied  with  the  strain  gages 


every  test  by  the  rapid  shunting  of  a  sequence  of  known  parallel  resistors 

with  the  response  recorded  on  an  oscilloscope  providing  the  scale  factor  for 

non 

the  strain 

IV. B.  EXPERIMENTAL  APPARATUS  AND  PROCEDURE 

The  experimental  apparatus  employed  for  the  present  comminution  studies 
of  rock  by  mechanical  impact  consists  essentially  of  a  pneumatic  gun  resting 
on  a  heavy  steel  table,  a  wooden  stand  on  which  the  sample  and  aluminum  back¬ 
up  cylinder  are  placed,  a  high-speed  camera  with  a  framing  rate  continuously 
adjustable  from  20,000  to  1  million  per  second,  associated  triggering  devices, 
strain  measuring  transducers  and  associated  recording  equipment.  A  schematic 
of  the  arrangement  is  presented  in  Figure  3,  and  Figure  4  shows  a  photograph 
of  the  set-up.  The  compressed  gas  gun  features  a  1/4  inch  inside  diameter 
steel  barrel  with  a  length  of  approximately  54  inches  and  utilizes  either 
compressed  house  air  with  an  upper  limit  of  180  psi  or  nitrogen  as  the  driving 
mechanism.  The  latter  gas  is  withdrawn  from  a  standard  bottle,  initially 
compressed  to  a  pressure  of  2250  psi,  through  a  control  panel  to  a  reservoir 
behind  the  breech  that  is  isolated  from  the  barrel  by  an  electromagnetic 
valve  with  an  upper  limit  of  1500  psi  capacity.  An  electrical  noise  suppres¬ 
sor  was  employed  in  series  with  the  valve  switch  to  prevent  the  appearance 
of  large  electrical  transients  due  to  valve  activation  on  the  oscillographic 
records.  The  muzzle  velocity  of  the  1/4  inch  diameter  steel  spheres  employed 
as  strikers  in  the  present  experiments  is  controlled  by  the  chamber  pressure 
that  is  adjusted  by  a  regulator  on  the  panel.  Firing  of  the  projectile 
inserted  in  the  gun  muzzle  oocurs  by  opening  the  electromagnetic  valve, 
permitting  discharge  of  the  gas  in  the  reservoir  and  acceleration  of  the 
projectile;  subsequently,  the  valve  is  closed  automatically  to  prevent  waste 
of  the  bottled  nitrogen. 

Two  horizontal  slots  cut  into  the  walls  of  the  barrel  near  the  muzzle 
end  permit  the  passage  of  two  light  beams  at  a  fixed  distance  6  inches  apart 
that  are  focused  on  two  photo-diodes  whose  output  is  connected  to  a  type 
3A1  dual  trace  amplifier  and  monitored  on  a  Tektronix  type  565  oscilloscope. 
Interruption  of  the  beams  by  the  passage  of  the  projectile  generate  flat- 
topped  transients  with  a  rise  time  of  about  10-20  microseconds  that  can  be 
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Fig.  3  -  Schematic  of  the  Experimental  Arrangement 


Fig.  4  -  Photograph  of  the 
Experimental  Arrangement 


employed  to  reliably  measure  the  projectile  velocity.  A  typical  photograph 
of  the  oscilloscope  record  of  this  device  is  shown  in  Figure  5.  A  helium- 
neon  laser  beam  was  directed  through  the  slot  nearest  the  breech  onto  a  photo¬ 
voltaic  cell  that  produces  a  signal  sufficient  to  trigger  delay  time  units 
incorporated  in  the  Tektronix  565  dual-beam  oscilloscopes. 

The  first  delay  timer  provides  for  the  discharge  of  a  70  pf,  4. 8  kV  con¬ 
denser  bank  into  a  General  Electric  FT  524  flash  tube  activated  through  a 
[12] 

thyratron  tube  .  This  event  simultaneously  serves  as  a  shutter  and  a  light 
source  for  the  operation  of  the  Beckman-Whitley  WB-2  camera  which  has  79 
independent  frames.  The  speed  of  the  rotating  prism  in  the  camera  is  controlled 
by  a  rheostat  and  is  displayed  on  an  electronic  revolution  counter.  A  para¬ 
bolic  reflector  collects  the  light  and  passes  it  through  a  sand-blasted  plate 
of  plexiglass  which  acts  as  a  diffuser  and  produces  a  nearly  uniform  field 
of  illumination  on  the  target.  The  second  delay  timer  acts  as  the  trigger  for 
the  remaining  five  oscilloscope  channels  so  as  to  capture  the  start  of  the 
strain  gage  signal. 

The  delay  setting  was  critical  due  to  the  extremely  short  time  interval, 

200  Ms,  during  which  the  initial  portion  of  the  strain  signal  could  be  viewed 
with  an  acceptable  resolution  relative  to  the  delay  time  interval  of  about 
1500  ps  for  the  camera  and  1600  ps  for  the  gages.  In  order  to  properly  deter¬ 
mine  the  magnitude  of  this  delay,  a  device  was  developed  consisting  of  two 
pieces  of  aluminum  foil  2  inches  square  which  were  placed,  one  on  each  side, 
over  a  1  inch  square  hole  cut  in  a  piece  of  1/16  inch  thick  plexiglass.  When 
this  foil  was  struck  by  the  projectile  in  a  preliminary  test  designed  to  set 
the  proper  delay,  the  two  pieces  of  foil  touched  and  completed  the  circuit  to 
produce  a  voltage  spike  from  a  suitable  generator. 

The  agar  plate  and  surrounding  plastic  tube  arrangement  previously 
employed  to  capture  the  ejecta  from  the  target  plate  was  replaced  by  a  system 
consisting  of  a  thin  polyethelyne  liner  placed  around  the  target  and  the  gun 
in  such  manner  as  to  completely  encase  the  region  initially  traversed  by  the 
ejecta.  This  liner  was  positioned  inside  a  rectangular  metallic  catcher  tank 
covered  on  the  interior  with  foam  sheets  that  served  as  a  cushion  to  reduce 
additional  fragmentation  of  ejecta  upon  contact.  The  polyethylene  sheets  were 
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Velocity  Measurement 
Sweep:  lOOpsec/div 

Delay  Measurement 
Sweep:  20psec/div 
Delay  Time:  1760  psec 


Velocity  Measurement 
Sweep:  lOOpsec/div 

Delay  Measurement 
Sweep:  20psec/div 
Delay  Time:  1780  psec 


Fig.  5  -  Typical  Photographic  Records  of 
the  Initial  Velocity  and  Delay  Measurements 
for  Run  D180 
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stretched  over  a  thin  metallic  frame  contained  within  the  tank.  This  arrange¬ 
ment  was  not  totally  effective  in  capturing  the  entire  mass  ejected  from  the 
target,  since  some  of  the  particles  perforated  the  container  lining  and 
embedded  themselves  in  the  foam  from  which  they  could  not  be  readily  recovered. 
The  tank  was  supplied  with  two  viewing  ports  that  permitted  photography  of 
the  comminution  event. 

The  collision  of  the  striker  with  the  target  produced  comminution  of 
the  rock  in  various  degrees,  accompanied  by  crater  formation,  penetration  or 
both,  depending  on  the  nature  of  the  target  material.  When  the  projectile  is 
not  embedded,  as  occurred  invariably  in  the  lower  velocity  tests  conducted 
more  recently,  the  striker  rebounds  with  diminished  velocity  after  a  contact 
period  of  the  order  of  microseconds,  simultaneously  with  any  ejecta  produced. 
Concurrently,  energy  is  transferred  to  the  interior  of  the  sample  in  the  form 
of  strain  waves  that  could  produce  damage  in  regions  remote  from  the  contact 
point. 

The  experimental  procedure  diverged  somewhat  depending  upon  whether  dio- 
rite  or  shale  targets  were  employed.  In  the  former  case,  the  composite  system 
including  the  aluminum  backup  plate  was  placed  on  a  contoured  wooden  stand 
inside  the  catcher  tank  with  the  impact  face  positioned  within  the  viewing 
area  permitted  by  the  ports.  The  shale  targets,  which  consisted  of  single 
potted  and  taped  disks,  were  supported  in  a  similar  fashion  except  that  they 
were  not  greased  to  the  back-up  plate,  but  clamped  to  it.  In  view  of  the 
lack  of  significant  ejecta  formation  for  this  target,  the  trajectory  of  the 
debris  from  this  type  of  specimen  was  not  determined  from  photographs  during 
the  second  year  of  operation.  Typical  high  speed  photographs  were  presented 
in  Reference  [8]. 

Thus,  the  tests  on  the  diorite  targets  involved  the  simultaneous  deter¬ 
mination  of  the  initial  striker  speed  from  time  interval  data,  the  final 
projectile  velocity  and  debris  kinematics  by  means  of  high-speed  photography, 
crater  topography  measurement  by  means  of  feeler  gages,  ejecta  size  distri¬ 
bution  by  electronic  scanning  methods,  elastic  wave  energy  in  the  target 
from  internal  strain  histories  and  crack  production  in  the  sample  by  post 
mortem  sectioning.  On  the  other  hand,  for  the  shale,  only  impact  velocity  of 
the  projectile  and  penetration  depths  were  determined. 
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For  the  photographic  measurements,  the  high  speed  camera  was  located  66- 
1/2  inches  from  the  impact  area  and  recorded  the  data  on  a  5  foot  length  of  35 
mm  Tri-X  film  which  was  transferred  into  the  stationary  film  holder  from  the 
storage  cassette.  The  data  obtained  in  this  fashion  includes  the  history  of 
the  striker  before  impact,  the  contact  of  the  bail  with  the  rock,  the  dis¬ 
placement  history  of  the  rebounding  projectile,  and  the  trajectory  of  the 
debris  generated  for  a  period  of  about  600  ys,  thus  involving  a  framing  rate 
of  132,000  per  second.  The  ejecta  consists  of  a  cloud  of  fine  particles, 
appearing  in  the  form  of  a  conical  shell,  as  well  as  of  some  larger  fragments. 
The  area  appearing  in  the  photographs,  between  the  sample  and  the  end  of  the 
gun  barrel,  lies  in  the  vertical  plane  vernal  to  the  target  and  is  about  4- 
1/4  inches  high  by  2-1/8  inches  wide. 

The  frontal  position  of  the  cloud  was  traced  along  the  extreme  rays  and 
the  velocity  of  the  fastest  particle  was  ascertained.  With  certain  assump¬ 
tions,  this  permits  the  evaluation  of  the  kinetic  energy  of  the  fragments  and 

the  rebound  energy  of  the  ball.  The  former  was  facilitated  with  the  aid  of 

[8] 

a  computer  program  that  requires  as  its  input  the  position-time  information 
contained  in  the  photographs.  Development  of  the  film  occurred  in  a  Microdol 
X  solution  diluted  in  the  ratio  of  1: 3  to  provide  minimum  graininess  in 
enlargements  of  the  photographs.  Such  blow-ups  were  produced  on  polycontrast 
paper  with  a  magnification  of  6.5  to  7,  and  all  required  information  was 
determined  from  these  enlargements. 

Concentric  circles  of  1/8,  1/4,  3/8,  1/2,  3/4  and  1  inch  diameter  were 
drawn  with  india  ink  about  the  contact  point  of  the  diorite  and  with  pencil 
on  many  of  the  shale  targets.  For  diorite,  the  volume  and  the  net  surface 
area  of  the  craters  were  determined  by  depth  measurements  using  feeler  gages 
over  a  10  x  10  mesh,  that  were  recorded  on  a  Tektronix  D13  dual  beam  oscillo¬ 
scope  with  AM  502  and  5A22N  differential  amplifiers.  These  data  served  as  the 

r  si 

input  for  a  computer  program1*  that  calculated  these  variables  and  also 
plotted  a  contour  map  which  provides  a  reasonable  three-dimensional  repre¬ 
sentation  of  the  crater. 

In  the  plurality  of  the  low-velocity  tests  using  a  1/4  inch  diameter 
steel  sphere  as  a  projectile,  the  crater  produced  in  the  shale  samples 
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consisted  of  a  hemispherical  indentation  with  short  (i.e.,  less  than  striker 
radius)  cylindrical  stems  emanating  to  the  surface.  In  some  cases,  a  shallow 
dish  was  additionally  formed  at  the  intersection  of  the  stem  and  the  target 
face.  A  few  tests  employing  larger  projectile  diameters  resulted  in  severe 
cracking  of  the  entire  sample  and  surrounding  material  and  were  thus  discon¬ 
tinued.  In  view  of  the  minor  role  of  the  dish  and  the  brittleness  of  the 
material,  only  maximum  penetration  measurements  were  executed  by  means  of 
feeler  gages  in  the  shale  targets. 

The  diorite  ejecta  trapped  in  the  catcher  tank  was  carefully  transferred 
to  a  petri  dish  and  weighed.  With  the  assumptions  that  the  density  of  the 
ejecta  is  identical  to  that  of  the  material  in  the  original  regions  of  the 
crater  and  that  no  compaction  takes  place  in  the  target,  a  mass  scale  factor 
was  computed  that  related  the  material  recaptured  to  the  total  ejected;  as 
explained  previously,  some  particles  escaped  entrapment  and  others  might 
have  been  lost  during  the  transfer  process.  Since  the  shale  exhibited  signifi¬ 
cant  permanent  deformation,  it  was  not  possible  to  utilize  this  procedure  for 
this  material.  The  shale  ejecta  was  also  collected  and  weighed,  but  its  size 
distribution  was  not  analyzed. 

2 

The  collected  diorite  particles  were  distributed  randomly  on  a  5961  ram 
glass  plate  that  was  then  subjected  to  mild  agitation.  The  plate  had  been 
coated  with  an  anti-static  spray  so  that  conglomeration  of  the  particles  and 
concentrations  near  the  perimeter  of  the  viewing  area  would  be  avoided.  A 
region  amounting  to  34.6  percent  of  the  plate  area  was  surveyed  by  an  elec¬ 
tronic  scanner,  a  Quant imer  720  unit  located  on  the  Davis  campus  of  the 
University  of  California.  This  device  is  capable  of  examining  specified 
rectangular  fields  containing  aggregates  of  particles  of  diverse  shapes  and 
sizes;  in  the  present  instance,  five  17.9  mm  x  22.26  mm  domains  were  analyzed 
within  the  plate  region  for  each  of  which  the  number  of  particles  contained 
in  ten  distinct  intervals  of  their  projected  area  were  enumerated.  In 

contrast,  in  the  first  year  of  operation,  several  hundred  areas,  each  of 
2 

8.15  ram  size  were  sampled  by  the  Quantimet.  This  information  was  scaled 
up  both  by  the  area  factor  consisting  of  the  ratio  of  total  plate  area  to 
that  electronically  scanned,  and  by  the  mass  factor  previously  defined. 
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Simple  particle  geometries  capable  of  being  described  by  a  single  dimension, 
selected  from  microscopic  examination  of  the  ejecta,  were  utilized  to  compute 
surface  areas  and  volumes  of  the  fragmented  crater  mass  based  upon  the  elec¬ 
tronically  measured  distributions. 

In  general,  the  sweep  rate  for  the  strain  gages  was  chosen  as  20  ps/cm 
(or,  in  one  case  10  us/cm),  a  compromise  between  adequate  resolution,  suffi¬ 
cient  length  of  the  trace  and  triggering  considerations.  The  amplification 
varied  from  20  to  100  mv/cm  so  as  to  permit  the  recording  of  two  traces  on 
a  single  photograph  without  interference  and,  further,  to  avoid  a  ratio  of 
pulse  height  to  width  larger  than  2:1.  Typical  data  from  these  transducers 
are  presented  in  Figure  6.  These  traces  can  be  employed  to  compute  the  wave 
energy  on  the  hypothesis  that  the  pulse  is  hemispherical  and  propagates  in 
a  homogeneous,  elastic  isotropic  medium  without  reflection  from  free  surfaces. 

The  present  report  also  addresses  itself  to  an  examination  of  the  internal 
damage  below  the  surface  of  the  crater,  a  topic  not  considered  in  the  previous 
report.  While  the  activity  of  the  last  year  has  shown  that  the  crack  network 
on  a  typical  diametral  section  can  be  displayed  for  measurement  purposes,  it 
proved  virtually  impossible  to  separate  the  cracks  produced  as  the  result  of 
the  impact  from  those  existing  previously.  Nevertheless,  such  an  examination 
was  executed  to  establish  bounds  on  the  new  surface  produced  as  the  result  of 
the  collision  that  consists  both  of  ejecta  and  crack  formation. 

Thus,  the  5.5  inch  diameter  target  plates  were  sectioned  through  the 
maximum  crater  diameter  using  a  0.08  inch  thick  water-cooled  diamond  saw; 
this  section  was  further  cut  to  a  width  of  2.5  inches  to  permit  machining.  The 
region  on  this  section  below  the  crater  was  then  lapped  using  progressively 
finer  grades  of  emery  paper  from  #320  to  4/0  grit.  The  samples  were 

then  placed  under  illumination  by  a  carbon  arc  on  a  Unitron  metallurgical 
microscope  with  a  magnification  of  13.7.  Polaroid  film  55  P/N  was  used  to 
record  overlapping  sections  of  the  specimens  and  a  composite  photograph  of  the 
damaged  area  was  produced.  All  visible  cracks  occurring  in  the  entire  central 
diametral  region  were  approximated  as  straight  lines  whose  actual  terminal 
coordinates  were  digitized  by  means  of  a  plotting  device  controlled  through  a 
program  placed  on  a  PDP-7  computer.  The  output  was  processed  by  means  of  a 


31 


r 


■ 


Fig.  6  -  Strain  Gage  Histories  as 
Measured  for  Run  D180;  Sweep  Rate: 
20  psec/div  for  all  gages 
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program  which  calculated  the  cracked  surface  area  based  on  the  hypothesis  of 
an  axially-symmetric  distribution  of  cracks;  a  plotting  device,  Calcomp, 
associated  with  the  CDC  6400  computer  located  on  the  Berkeley  Campus  of  the 
University  of  California,  provided  an  automatic  graphical  representation 
of  this  crack  network. 

The  polished  crack  samples,  an  example  of  which  appears  in  Fig.  7a,  were 
then  dyed  using  Spotcheck  dye  penetrant  manufactured  by  Magnaflux  Corporation, 
Chicago  after  cleaning  of  the  section  with  a  solvent;  the  penetrant  was  applied 
as  a  spray  and  removed  with  a  cloth  after  30  minutes  in  the  case  of  the  fine 
cracks  present.  The  developer  was  then  sprayed  on  the  surface  and  allowed  to 
act  for  a  period  ranging  from  5  to  15  minutes  before  inspection  of  the  surface; 
the  cracks  and  pores  are  brought  out  as  a  series  of  red  lines  and  dots.  This 
pattern  was  photographed  with  55  P/N  Polaroid  film  and  enlarged  to  the  same 
magnification  as  the  composite  picture  of  the  damaged  area  in  the  undyed 
specimen,  permitting  an  assessment  of  the  regions  and  amount  of  damage.  In 
a  similar  manner,  an  area  of  the  section  remote  from  the  crater  (at  least 
one  crater  diameter  away)  was  surveyed  on  the  hypothesis  that  this  region 
was  representative  of  the  virgin  state  of  the  rock.  The  difference  between 
these  networks  was  attributed  to  the  impact  process. 

Several  quasi-static  tests  were  conducted  on  the  shale  with  an  Instron 
machine  using  various  rates  of  feed.  A  4  inch  long,  1/4  inch  diameter  section 
of  drill  rod  with  a  hemispherical  tip  to  simulate  the  shape  of  the  striking 
sphere  was  pressed  into  the  sample  to  obtain  a  force-indentation  relation  and 
to  determine  if  strain-rate  effects  existed  in  the  quasi-static  range.  The 
indentation  rate  was  maintained  constant  and  the  force  history  was  recorded  by 
means  of  a  load  cell.  These  data  permit  the  calculation  of  the  energy  of 
penetration  as  a  function  of  depth. 

IV. C.  RESULTS  AND  DISCUSSION 
1.  Diorite 

a.  Velocity  Measurements 

For  the  most  part,  the  initial  velocity  of  the  projectile  as 
determined  by  the  photodiode  measurement  was  higher  than  that  calculated 
from  the  high  speed  camera  data;  hcwever  for  two  different  pressures, 
the  camera  data  yielded  the  same  velocity. 
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Since  some  of  the  runs  yielded  an  insufficient  number  of  frames  for 
reliable  initial  velocity  measurement  from  the  photographs,  the  original 
projectile  speed  was  uniformly  obtained  from  the  light  beam  interruption 
technique.  The  camera  data  was  considered  to  be  reliable  for  the  deter¬ 
mination  of  the  ejecta  and  projectile  rebound  velocities.  The  velocity 
of  the  leading  edge  of  the  cloud  formed  by  the  ejecta  mass  is  given  in 
Table  5. 

b.  Ejecta 

A  typical  high-speed  photograph  of  the  ejecta  cloud  is  shown  in 
Figure  8.  The  kinetic  energy  of  these  particles  was  calculated  on  the 
assumptions  that  elemental  rings  of  constant  mass,  dm,  are  emitted  during 
impact  in  the  form  of  a  right  circular  conical  shell  whose  axis  is  taken 
along  the  x-direction,  and  that  the  velocity  of  each  ring,  v,  varies 
linearly  from  a  maximum  measured  value  vm  at  the  cloud  front  to  zero  at 
the  impact  face,  as  shown  in  Figure  9  which  represents  the  instant  of 
final  ejecta  formation  when  the  altitude  of  the  shell  is  L.  These 
assumptions  are  considered  to  be  reasonable  in  view  of  previous  analyses 

r  s  3 

of  the  spatial  particle  distribution1  J  and  inspection  of  the  present 
photographs  that  indicate  the  cone-like  shape  of  Lhe  cloud,  and  that, 
further,  successively  emitted  particles  move  at  a  consecutively  slower 
rate.  In  fact,  some  ejecta  were  found  to  be  resting  in  the  crater  upon 
inspection  after  impact. 

The  energy  associated  with  the  element  of  mass  dm  is  given  by 
dEej  =  (1/2) v2  dm  (2) 

The  symmetry  of  the  cloud  and  the  velocity  distribution  assumption  leads 
to  the  relations 

v  =  vm(x/L)  and  dm  =  (M/L)dx  (3) 

where  M  is  the  total  ejecta  mass,  determined  on  the  basis  of  the  measured 

crater  volume.  Upon  substitution  and  integration  of  Equation  2,  the 

total  kinetic  energy  of  the  ejecta,  E  .  becomes 

ej 

Eej  = 
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TABLE  5.  EJECTA  CLOUD  CHARACTERISTICS 
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graphs;  9^  determined  from 
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Fig.  9  -  Assumed  Ejecta  Configuration  and  Velocity 
Distribution 


While  this  appears  to  be  a  reasonable  estimate,  it  is  considered  to  be 
a  lower  bound  value  in  view  of  the  prescribed  velocity  distribution. 
However,  this  hypothesis  is  much  more  reasonable  than  the  use  of  a  con¬ 
stant  value  of  the  velocity  v^  for  each  particle  of  the  cloud,  as  em¬ 
ployed  in  Reference  8,  although  the  truth  is  probably  somewhere  in 
between,  closer  to  the  current  assumption.  It  is  not  possible,  at  the 
present  time,  to  measure  the  velocity  variation  within  the  particles  of 
the  densely  packed  cloud.  Furthermore,  the  kinetic  energy  of  the  ejecta 
is  a  small  percentage  of  the  initial  energy,  ranging  from  7  to  15  per¬ 
cent  in  the  present  tests. 

In  the  investigation  of  the  first  year,  the  kinetic  energy  of 
all  fragments  of  sufficient  size  to  be  individually  identified  in  the 
high-speed  photographs  was  computed  separately.  This  procedure  was 
not  followed  in  the  present  calculations  because  (a)  the  quantity  of 
energy  involved  was  small  relative  to  the  cloud  energy,  (b)  significant 
inaccuracies  in  the  magnitude  calculated  for  the  energy  of  these  special 
fragments  resulted  from  the  lack  of  knowledge  of  the  inclination  of 
the  trajectory  with  the  normal  to  the  specimen  face,  and  (c)  in  some 
cases,  it  was  difficult  to  identify  a  particular  fragment  in  a  sequence 
of  frames.  Thus,  all  ejecta  was  represented  by  a  single  distribution 
function. 

The  results  of  the  investigation  during  the  first  year  had 
provided  a  number  of  curves  of  particle  size  as  a  function  of  spatial 
distribution,  that  showed  a  Gaussian  trend  for  the  diorite;  no  further 
tests  of  this  nature  were  considered  necessary.  On  the  other  hand,  the 
particle  size  distribution  in  each  of  the  tests  was  electronically 
determined  from  the  sampled  regions.  A  typical  photograph  of  the  par¬ 
ticles  caught  on  the  agar  plate  is  reproduced  in  Figure  10  that  shows 
the  projected  particle  shapes  of  the  ejecta;  a  reasonable  approximation 
of  these  geometries  appear  to  be  circles  and  ellipses.  The  data  for 
the  number  of  particles  vs.  projected  surface  area  obtained  from  the 
Quantimet  samplings  were  fed  into  a  computer  program  that  scaled  the 
mass  of  the  particles  collected  to  that  of  the  crater  using  the  same 
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Fig.  10  -  Typical  Quantimet  720  Unit 
Photographs;  2.5  mm  x  3.2  mm  Field  of 
View 


4l 


■ 


particle  size  distribution  as  that  experimentally  determined.  The  output 
consisted  of  the  number  of  particles,  fraction  of  total  particles,  cumu¬ 
lative  surface  area  and  cumulative  volume  of  the  particles  for  each  size 
range,  based  on  two  assumed  simple  particle  geometries  described  by  a 
single  dimensional  parameter.  The  shapes  hypothesized  were  (a)  an  oblate 
spheroid  with  a  series  of  ratios  of  minor  to  major  axis  lengths  ranging 
from  0.1  to  1.0  in  0.1  increments,  and  (b)  disks  with  a  ratio  of  height 
to  radius  of  either  0.3  or  5.0.  These  geometries  conformed  quite  well 
to  the  observed  projected  shapes  of  Figure  10.  A  typical  single  test 
result  for  each  of  the  three  assumed  shapes  consisting  of  the  two  disks 
and  a  sphere  is  presented  in  Table  6.  A  histogram  of  the  fraction  of 
total  particles  derived  by  scaling  the  distribution  obtained  from  the 
Quantimet  readings  for  the  particles  actually  collected  to  that  corre¬ 
sponding  to  the  entire  crater  volume  as  a  function  of  square  picture 
points  is  shown  in  Figure  11  for  an  assumed  spherical  particle  shape. 
Square  picture  points  represent  the  unit  of  viewing  area  in  the  Quantimet 

sampling  that,  for  the  present  magnification,  represent  a  domain  8.24  x 
-4  2 

10  mm  .  The  cumulative  information  derived  from  Figure  11  has  been 
replotted  on  a  probabilistic  scale,  as  indicated  in  Figure  12,  yielding 
a  good  approximation  to  a  normal  fit.  The  same  conclusion  was  drawn 
from  the  results  of  Reference  8  where  the  particles  were  collected  on  a 
gelatinous  catcher  plate  and  the  size  and  spatial  distributions  were  both 
evaluated.  The  logarithmic  mean, Lny,  and  standard  deviation, a  ,  associated 
with  a  Gaussian  description  of  the  size  distribution  has  been  summarized 
in  Table  7  as  a  function  of  input  energy  for  three  specific  particle 
shapes,  and  the  logarithmic  mean  has  also  been  plotted  in  Figure  13. 

The  discontinuities  in  this  diagram  are  the  result  of  the  large  varia¬ 
tions  found  between  samples,  although  the  average  can  be  readily  drawn. 
From  the  data  collected,  there  appears  to  be  a  trend  towards  decreasing 
particle  size  with  increasing  energy. 

As  may  be  noted  from  the  summary  for  the  diorite  runs,  Table 
8,  the  ratio  of  recovered  particle  mass  to  crater  mass  ranged  from  about 
40  to  80  percent.  It  was  assumed  in  the  previous  calculations  that  this 
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TABI£  6a 


Gt-liWL  T K iLAt.'  PKLHLKULS  Of  iHt  EJLCTA  -  UIORITE  1«8 

PARTICLE  SHApt  ASSUMED  10  BE  A  C1SK 
HtlbMT  s  S.B  RADIUS 


AREA-CROSS 
SECI 10RAL 
PIC. PIS 

ARLA-bWE 
PAM  IClL 
SO  .  MM 

VUL.-UNL 
PAR  1  IDLE 
SU.«M 

NLMBt  R  'OF 
PARTICLES 

FRACTION 

particles 

area-all 

particles 

SQ.MM 

VOL. -ALL 

particles 

CU.MM 

l.BBdt+kJl 

3.bb»L-l<32 

1.4blt  +  k)l 

i.4«er-t»2 

5.34SE-01 

9.122L-03 

7.3i9E-b2 

i.  /fab£.-HA 

i.niLt*2 

1.630E-01 

1.245E+01 

3 , 004L-01 

4.k5k)k)L+.'l 

1.4b4L-01 

4.99bE-k>A 

1.744L+02 

1 .672E-0 1 

2.554E+U1 

8 , 71bE-l41 

d.MkJtft  +  01 

2«92bt-ul 

3  •  3fo7L  +  kJ2 

3.17HE-01 

9.683E+01 

4.R7HL+00 

l.bltfBEW 

tj.abbL-kJl 

3.99/ L-Mi! 

1.492E  +  KI2 

1.43BE-B1 

6.73bE+01 

5.9bbL+00 

3.2k!0L+i)2 

1  •  1  /IL  +  t  i) 

l.Hfa5L  +  t)2 

1.U4BE-01 

1.271E+02 

1.22/L+01 

.  fa.4k)tL-*-t2 

2.3fak)L  +  uk) 

3.23bF-Rl 

b.9b0t+ul 

6.690E-B2 

1 . 64BE+02 

2 . 2SBL+01 

1.3k)kJL  +  «i3 

4  ,  bbDL  +  dtf 

a.9dbE-«i 

l.AblL+Ul 

1.40BE -0^ 

6.814E+01 

1.315L+B1 

<;,SkJ*JE+*'3 

9.  149t  +  dt) 

2.4691 +0« 

b.“69t  +  k)y 

6.200E-03 

S.91faC*01 

1.597L+01 

l.bbPE+Dl 

b.9tt2E+00 

4.909E+HH 

4.7PBE-03 

a.97bE+01 

3.424E+01 

1  ,HD0E+ki4 

TOTAL  (PaRT«  FKAC I  * 

AREA i  VUL» 

l,K43t  +  li)3 

i.kJPHE  +  0‘1 

7 , 316E+02 

l,lk)0t  +  »2 
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TABLE  6b 

GEOMETRICAL  PKUHEKTitS  OF  IrL  EJECTA  -  UIORITE  lad 

PARTICLE  ShAJ»E  ASSUMED  1 0  BE  A  CISK 
HElbHT  =  3/10  RADIUS 


AREA-CROSS  AREA-ORE  VUL.-UNt 
SECTIONAL  RAH  I  IDLE  par  hole 
HIC.ETS  SQ.MM  SCI. MM 


NUMBER  OF  FRACTION  AREA-ALL  VOL.-ALL 
PAHIIClES  particles  particles  particles 

SQ.MM  CU.MM 


_i.iUii>L+0l 
2 .  0  00E  +  01 
4.0B0E+01 
_e.j0HRL+iJi 
1.600L+02 
3 .200L+il)2 
_ii»HIK0L+iJ2 
1.300E+03 
2.b00L+d3 
baWHWEtL3. 


3.029E-02  2.128E-04  4.287E+01  1.400E- 
£i • 0b7E-U2  fa.Ml9E-MR  4.991E+02  1.63PE- 
1.211E-01  l./»2E-k)A  S.  120E+02  1.672E* 
2.423L-01  4.aibE“0A  9.706E+02  3.170E- 
4. 84bt-.nl  1.3fa2E*02  4.379E  +  E2  1.43BE' 
9,£>92L-0l  3.ab^E-0*  3.164E+02  1.040E 
1.9b3E*00  1.102E-01  2.H48E  +  02  6.69BE 
3,  atilt +  00  3,0bAE-0i  4.2b7t+0l  1.40HE 
7.572L+00  a.*+i2E”0l  l.a98E+0l  6.2R0E 
1.514E*0i  2.3/9E+00  l.439E+0l  4.700E 


02  1.298E+00 
■01  3.025E+U1 
01  6.202E+B1 
•01  2.352E+H2 
.01  2.122E+02 
.01  3.H8faE+02 
-02  4.001E+02 
-02  1.655E+02 
-03  1.437E+02 
-03  2 . 179E+02 


9.122L-03 

3.00HE-01 

8.71SL-01 

4.674E400 

5.963E400 

1.227L+01 

2.2b8E+0l 

1.31AE+01 

1.597L+01 

3.424E+01 


1.000E+U4 

TOTAL  (PART*  FRACm  AREA*  vuli 


3.0&2C+03  1.000E+00  1.777E+03  1.100E+02 
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TABI£  6c 


GtOHLTKiLAL  PKUPLKTILS  Ul-  1  Ht-  EJLCTA  -  DIORITE  180 
MINUH  AXl^/MAJOR  AXIS  =  1.000  -  -  - 


AKtA-CkUoS 
SECTIONAL 
PIC. HIS 

AKL  A”ui'lC 
PAAl  lCLfc. 
..  Su.MM. 

VUL.-UNL 
PAH  1 lLLk 
SQ.MM 

NUMBER  OF 
HAKTICLtS 

rRACTlCIM 

PAHTIClES 

area-all 

particles 

SQ.MM 

VOL. -ALL 

part iclls 

CU.MM 

l.K)00fc.^~l 

4,6S>9E-0*L  9,Ht.8L-0‘+ 

9.6H5E+00 

1.400E-02 

H.494E-01 

9.122E-03 

2.k)0k)L  +  iJl 

9.319L-02 

2.b/bE-03 

1.123E+02 

1.630E-B1 

1.047E+01 

3.004L-01 

4.00k>L+*ll 

1 ,8b4L-kJl 

7.bb6L-03 

1.152E+B2 

1.672E-01 

2.147E+B1 

8.715L-01 

8.k)00L  +  kil 

3.728L-01 

2.1H0E-02 

2  •  lb4E+02 

3.170E-01 

8.141E+01 

4.674L+00 

l.bk>0L  +  u2 

7.4b5t-Hl 

6.B53E-02 

9.852E+01 

1.430E-B1 

7.345E+01 

5.963L+00 

3.200t+i:2 

1.491E+00 

1.  /I2E-01 

7.1fc5E*0l 

1.040E-01 

1.068E+02 

1.227L+01 

6.AkJ0t  +  H2 

3,k)k)bt  +  k>0 

4,8y>3E-0l 

4.b09L+01 

6.690E-B2 

1.385E+H2 

2 .2SBE  +  01 

1.A00L+M3 

b.94tfE  +  k1k) 

l.bbIE+00 

9.645E+0B 

1 .40PE-02 

5.729E+01 

1.313L+01 

2.bk)0E  +  03 

1  •  lbbL-f  0 1 

6, /38E+H0 

4.271E+M0 

6.200E-03 

4.97&E+01 

1.597L+B1 

b.000L+03 

2.3Ak)t*«l 

1.0b /E+01 

3.23BE+00 

4.700E-B3 

7.544E+01 

3.424t40l 

1 ,  kJkiOL-*  04 

TOTAL  ( PART  «  FIvACI, 

A  R  E  A  t  VOL) 

6.tib9L  +  ki2 

1.000E+0H 

6.151E+B2 

1.100L+02 
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Ln/t,  LN  MEAN  PARTICLE  SIZE 

Fig.  13  -  Energy  Input  versus  Mean  Particle  Size  for 
Three  Particle  Shapes 


TABLE  8.  DIORITE  TEST  RESULTS 
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Photodiode  velocity  measurement 


TABLE  8  (cont'd).  DIORITE  TEST  RESULTS 
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Based  on  a  spherical  particle  shape 


TABLE  8  (cont’d).  DIORITE  TEST  RESULTS 


non-recovered  mass  followed  the  same  size  distribution  as  that  collected 
in  the  catcher  unit.  This  assumption  is  probably  in  error,  since  the 
fragments  most  likely  to  have  escaped  the  measurement  process  were  those 
at  the  extreme  end  of  the  distribution  spectrum.  The  larger  particles 
had  the  best  chance  to  perforate  the  liner,  while  the  smallest  ones, 
amounting  to  dufit,  were  most  likely  lost  in  the  transfer  operations. 

Thus,  an  appraisal  of  the  effect  of  these  losses  on  the  actual  particle 
size  distribution  would  indicate  that  the  mean  value  was  probably  not 
strongly  affected,  but  that  the  deviation  might  be  considerably  larger 
and  the  distribution  thus  widened.  The  total  ejecta  surface  area  is 
most  likely  not  changed  a  great  deal  by  this  increment  in  the  distribu¬ 
tion. 

c.  Wave  Energy 

In  a  previous  parallel  investigation,  Reference  [13],  an  esti¬ 
mate  was  given  of  the  amount  of  energy  in  the  hemispherical  wave  pro¬ 
duced  in  a  block  of  homogeneous,  isotropic  elastic  material  due  to  a 
prescribed  time-dependent  load  applied  over  a  small  portion  of  the 
surface  of  the  material.  The  values  obtained  were  based  on  calculations 
of  an  exponentially  decaying  pressure  input  applied  to  the  surface  of 

rx4] 

a  small  spherical  cavity  on  the  inside  of  an  infinite  medium  and  a 
prescribed  ratio  of  tangential  to  radial  strain,  assumed  to  be  1/2. 
Surface  strain  gages  provided  the  information  for  the  value  of  peak 
radial  strain  that  was  taken  to  be  the  maximum  pressure  in  the  solution 
of  the  boundary  value  problem.  The  result  indicated  that  the  energy 
fraction  consumed  in  the  elastic  wave  amounted  to  from  about  5  to  15 
percent  of  the  initial  energy  of  the  projectile.  The  present  effort 
involves  the  measurement  of  both  radial  and  tangential  strains  in  the 
interior  of  the  specimen  and  a  more  accurate  analytical  representation 
of  the  physical  process  in  an  attempt  to  obtain  a  better  evaluation  of 
the  magnitude  of  this  energy  component. 

The  wave  energy  was  calculated  on  the  basis  of  the  identical 
basic  assumption  as  above  of  a  perfectly  spherical  elastic  transient 
emanating  from  the  impact  point  through  one-half  of  an  infinite  space 
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of  density  p.  At  any  distance  r  from  the  contact  location,  the  work  done 

W  by  the  pressure  -  cr  on  the  internal  hemispherical  cavity  is  given  by 

rr  T 

W  ■  -  2  t  r  f  o  <t>u  <t>  dt  (5) 

o  rr  v 

where  u  is  the  radial  displacement,  t  is  time  measured  from  the  start  of 

the  pulse  at  position  r,  t  is  the  pulse  duration,  and  a  superposed  dot 

represents  differentiation  with  respect  to  time.  The  radial  stress  o 

is  given  from  Hooke's  law  as 

orr  =  (A  +  2G)err  +  A(r00  +  e^)  (6) 

where  the  principal  strains  in  directions  r,  0  and  are  given  by  z  , 

/  rr 

e00’  e^>  an(*  *  and  ®  are  ^ame  constants.  For  a  symmetric  wave, 

e00  =  an^  hence 


a  <t> 
rr 


(A  +  2G)err  <t>  +  2Ae00<t> 


(7) 


The  particle  velocity  u  is  then  given  by 


“  CD  (err  +  eee) 


where  c  *  dilatational  velocity 


■  JA+2G 

»  o 


Substitution  of  Equations  (7)  and  (8)  in  Equation  (5)  yields 


W-2ttc 


A 


T  T 

(A+2G)  f  z  2  <t>dt+(3A+2G)  /e 


rr  <t>Eee<t>dt+2A  {c6o‘ 


<t>dt 


(8) 


o  rr  '  o 

Since  the  radial  and  tangential  strains  were  measured  at  different  radii 
* 

r  and  r  ,  respectively,  it  was  necessary  to  obtain  their  equivalent 
value  at  the  same  position.  From  the  strain-displacement  equations,  this 
equivalence  is  given  by 


(9) 


'00 


(r*/r)2e* 


00 


(10) 


Substitution  of  this  relation  in  Equation  (9)  yields  the  wave  energy 

WEw-2TrcDr2  [(A+2G)  fzj  dt+(3A+2G)  (r*/r) 2  c^dt+2  A(r*/r)4  /ej02dt]  (11) 

in  terms  of  the  measured  strains  z  and  e*  whose  location  was  given 

rr  90  ° 

in  Table  3  and  whose  history  is  exemplified  in  Figure  5. 

A  computer  program  was  written  to  evaluate  the  integrals  of  Equation  (11) 
using  a  trapezoidal  rule. 

A  sample  problem  was  constructed  to  test  the  accuracy  of  the 
program  using  assumed  strain  histories  consistent  with  a  spherical  wave, 
given  by 


54 


0  £  t  £  T 

otherwise 


(12) 


'80 


e 

rr 


(uQ/r)  sin  — 
0 


TTt 

T 


+  (irr/CjjTjcos 


0  £  t  £  t 
otherwise 


Utilizing  test  data  for  u0  and  t  yielded  computed  values  of  Ew  that 
were  within  0.1  percent  of  the  result  obtained  by  exact  integration. 

For  two  of  the  tests,  a  malfunction  occurred  in  the  circuitry 
of  tangential  gage  T2.  This  demanded  a  further  investigation  to  deter¬ 
mine  whether  the  wave  energy  could  be  calculated  from  the  history  of 
only  one  radial  gage  with  the  necessity  of  additional  assumptions. 

For  an  outward  propagating  spherical  wave,  starting  at  time  t  =■  0, 


u  ■  b  -  cd>]a 

and  since  the  strain-displacement  equations  are  given  by 


3u 

Err  3r  * 
it  follows  that 


'00 


_u 

r 


(13) 


(14) 


Err<t*  r>=  "  72  F<t  -  *f->  - 


_r 

c„ 


1 

rcr 


F'<t — —  >and  e00<t,  F<t-  ^->  (15) 


D  D  '"D  vv  r ‘ 

where  the  prime  denotes  the  derivative  with  respect  to  the  argument. 
Upon  multiplication  of  Eqi 


b  •  b]  /r‘ 


'00 


relations,  it  follows  that 


uation  (15)  by  -  Cp/r,  use  of  the  identity 


F'<t 


-v] 


/rz,  and  the  strain-displacement 


cee  <l;>  +  T  Eee  <t;> 


—  e  <t> 
r  rr 


(16) 


Integrating  this  ordinary  differential  equation  with  the  initial  condi- 

0  for 


tion  e00<t>  *  0  for  t  =  0  (this  is  equivalent  to  e00<t;.  r> 
t  £  r/cD)  yields 


1 99 


CD  'CDt/r  ,C  _  ,  V/E 
—  e  /  e  <t>  e  dx 

r  ;  o  rr 


0  <  T  <  t 


(17) 


As  a  check  on  the  assumptions  inherent  in  the  derivation  of  Equation 
(17),  the  energies  were  calculated  for  all  cases  where  both  radial  and 
tangential  strain  histories  were  available,  and  the  results  computed 
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by  means  of  Equation  (17)  were  compared  with  the  data.  Reasonable  corre¬ 
spondence  was  obtained  as  indicated  in  Table  9  which  exemplifies  the  wave 
energy  calculations  for  Run  D180. 

The  wave  energy  was  computed  for  each  test  using  gage  R1  and  the 
shifted  value  of  gage  T2  as  the  principal  stations  for  the  radial  and  tan¬ 
gential  measurement,  respectively.  Since  the  tangential  gage  failed  in 
runs  D120  and  D240,  Equation  (17)  served  as  the  replacement  for  these  data. 

Figure  14,  which  shows  the  actual  and  calculated  values  of 
for  Run  D180  indicates  the  nature  of  the  discrepancies  in  the  amplitude, 
rise  time  and  duration  of  the  two  histories.  Although  these  differences 
are  not  excessive  and,  as  a  result,  it  is  not  expected  that  significant 
changes  in  the  wave  energy  would  result  from  these  variations,  they  do 
indicate  that  the  basic  assumptions  of  a  spherically-symmetric  wave  in  a 
homogeneous,  isotropic  elastic  medium  are  not  satisfied  exactly  in  the  pre¬ 
sent  instance,  as  known  a  priori;  however,  these  hypotheses  do  provide  re¬ 
sults  of  sufficient  accuracy  for  the  present  purpose,  particularly  since 
the  wave  energy  represents  a  relatively  small  portion  of  the  impact  energy. 

d.  Crater 

As  in  the  earlier  work,  crater  topography  measurements  were  exe¬ 
cuted  by  means  of  feeler  gages;  a  typical  contour  map  is  presented  in 
Figure  15  and  the  volume  and  surface  area  of  the  craters  are  summarized  in 
Table  8.  The  crater  contours  ordinarily  exhibited  at  least  one  axis  of 
symmetry  and  a  relatively  smooth  interior  shape.  The  crater  depth  and  the 

maximum  crater  diameter  were  approximately  one-half  and  one  projectile 

r8] 

diameter,  respectively,  as  was  also  observed  previously 

e.  Cracks 

Although  inhoraogeneities,  previously  existing  cracks,  stress 
concentrations  and  stored  internal  energy  might  alter  the  axial  symmetry 
of  the  response  of  the  rock  sample  to  the  central  impact  by  the  steel 
sphere,  the  influence  of  these  factors  cannot  be  quantitatively  evaluated 
in  the  present  or  any  foreseeable  investigation.  In  consequence,  as  also 
somewhat  reinforced  by  visual  inspection  on  a  diametral  section  of  the 
damage  pattern  in  the  specimens,  the  extent  of  the  crack  system  produced 
by  the  impact  had  to  be  assumed  and  for  the  sake  of  simplicity 
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TABI£  9a  (cont'd) 
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was  taken  to  be  symmetric  about  the  impact  axis,  f  ming  surfaces  of 
revolution.  A  computer  program,  whose  listing  is  given  in  Appendix  A, 
was  developed  that  calculated  the  area  associated  with  these  new  sur¬ 
faces,  the  length  and  angular  orientation  of  each  crack  and  their  total 
number  based  on  an  input  that  stipulated  the  end  points  of  the  cracks 
with  respect  to  a  coordinate  system  centered  at  the  point  of  impact;  it 
was  assumed  that  the  cracks  on  the  diametral  surface  were  straight  lines. 
These  inputs  were  also  employed  to  generate  a  computer  plot  of  the  dia¬ 
metral  section,  such  as  shown  in  a  five-fold  enlargement  in  Figure  16; 
the  corresponding  13.7  times  photographic  enlargement  of  the  actual 
specimen  is  presented  in  Figure  7.  The  scheme  of  the  machined  sections 
is  also  included  in  this  figure. 

The  crack  surface  so  computed  cannot  distinguish  between  that 
produced  by  the  impact  and  that  preexisting.  In  order  to  attempt  to 
account  for  the  initially  present  crack  surfaces,  the  computed  areas 
were  adjusted  on  the  basis  of  a  crack  network  sampled  in  two  selected 
specimens  D120  and  D220  at  regions  A  or  C  which  were  so  remote  from 
the  impact  point  that  they  were  believed  not  to  have  been  affected  in 
any  way  by  the  contact  phenomenon  and  were  thus  representative  of  the 
virgin  state.  The  crack  surface  per  unit  volume  believed  to  exist 
originally  throughout  the  specimen  was  computed  both  on  the  basis  that 
the  axis  of  symmetry  passed  through  the  center  of  the  region  examined, 
and  on  the  basis  of  an  isotropic  crack  pattern  using  the  average  value 
of  the  cracks  for  the  two  specimens  D120  and  D220.  The  net  crack  surface 
due  to  the  impact  was  taken  as  the  difference  between  these  values;  it 
is  shown  in  Table  10  for  both  methods  of  calculation.  It  is  clear,  in 
view  of  the  negative  crack  areas  calculated  in  this  fashion,  that  the 
procedure  adopted  is  highly  questionable,  even  on  an  average  basis. 
Apparently,  individual  variations  in  the  initial  crack  configurations 
grossly  overshadow  that  generated  by  the  impact  as  was  also  confirmed 
by  visual  inspection  of  all  samples  tested,  and  average  values  simply 
cannot  be  applied  to  the  present  situation. 
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Fig,  16  -  Typical  CALCGMP  Plot  Approximating 
the  Cross -sectional  Crack  Distribution 
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Another  attempt  was  made  to  overcome  these  individual  variations. 
The  entire  central  diametral  section  of  the  specimens,  domain  B,  was 
covered  with  the  dye  penetrant  and  the  developer  was  applied.  This 
yielded  disjointed  regions  of  visible  crack  damage,  as  shown  in  Figure 
17.  The  crack  network  in  a  subregion  of  the  damaged  domain,  a  distance 
of  at  least  one  diameter  away  from  the  crater,  was  examined  and  the  ratio 
of  the  dyed  (cracked)  to  the  total  surface  area  within  that  domain  of 
scrutiny  was  deLermined.  This  area  factor  was  now  applied  to  adjust 
the  initially  assumed  preexisting  crack  distribution,  determined  on  the 
basis  of  the  average  of  samples  D120  and  D220,  to  conform  to  the  cracks 
present  in  the  actual  cross  section  of  each  test  specimen  by  multiplying 
this  average  crack  surface  area  by  the  quotient  of  the  ratios  of  dyed 
to  total  areas  of  the  actual  and  the  two  selected  specimens.  The  result 
was  considered  to  be  the  estimate  of  the  preexisting  damage  in  each  of 
the  test  specimens,  and  is  also  shown  in  Table  10.  A  net  crack  surface 
generated  by  the  impact  based  on  this  preexisting  crack  computation  is 
also  included.  Tinally,  a  calculation  is  summarized  in  the  table  which 
bases  the  preexisting  crack  surface  area  on  the  total  crack  length  per 
surface  area  found  in  the  region  under  examinarion;  this  quantity  can 
readily  be  related  to  crack  surface  per  unit  volume.  As  may  be  seen, 
all  of  these  approaches  yielded  negative  crack  surface  estimates  pro¬ 
duced  by  the  impacts  for  at  least  some  of  the  specimens,  a  situation 
that  is  considered  to  be  phenomenologically  unrealistic.  This  situation 
exists  even  when  the  procedure  is  applied  to  the  same  central  section 
of  the  selected  samples  whose  crack  area  several  diameters  away  from 
the  impact  region  was  employed  as  the  standard  for  the  fabric  of  the  vir¬ 
gin  specimen. 

Since  it  is  not  possible  to  determine  the  extent  of  preexistirg 
cracks  prior  to  the  test  without  destroying  the  specimen,  it  is  either 
necessary  to  utilize  statistical  results  that  are  physically  acceptable, 
even  when  some  of  the  ingredient  information  violates  a  reasonable  inter¬ 
pretation  of  the  phenomenon,  or  else  be  satisfied  with  bounds  on  the  new 
surface  generated  that  ignores  the  crack  sampling  outside  the  impact  aone. 
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TABLE  10.  CRACK  SURFACE  AREA  BASED  ON  VARIOUS  ASSUMPTIONS 
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Adjusted  areas  for  all  runs  involve  the  average  of  D120  A  and  D  220  C  areas 


TABLE  10  (cont'd).  CRACK  SURFACE  AREA  BASED  ON  VARIOUS  ASSUMPTIONS 
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Remote  measurement  scaled  by  crack  length  ratio 


The  average  of  all  positive  and  negative  generated  crack  surfaces  based 

on  the  dying  procedure  and  application  of  the  area  factors  was  evaluated 
2 

as  1490  mm  for  the  impact  of  a  1/4  inch  diameter  steel  sphere  in  the 
velocity  range  from  169  to  246  m/sec;  individual  variations  within  this 
span  were  completely  obscured  by  the  differences  within  the  samples 
supplied. 

In  view  of  the  above  considerations,  the  energy  required  to 
produce  unit  new  surface  will  be  computed  on  the  following  bases: 

(a)  All  cracks  located  in  the  central  diametral  region  that  extended 

virtually  to  the  edge  of  segment  B  of  the  specimens  were  considered  to 

be  produced  by  the  impact,  (b)  A  net  crack  surface  based  on  the  average 

2 

new  surface  considered  to  be  produced,  amounting  to  1490  mm  was  incor¬ 
porated  in  addition  to  the  ejecta  and  net  crater  surfaces,  and  (c)  All 
observed  cracks  were  assumed  to  be  preexisting.  This  information  is 
compiled  in  Table  11,  and  is  the  energy  of  comminution  required  to 
produce  new  surface  area  based  on  the  hypotheses  indicated.  Obviously, 
the  scheme  adopted  here  is  such  that  the  magnitude  of  this  quantity 
increases  from  assumption  (a)  to  assumption  (c).  For  the  range  inves¬ 
tigated,  the  energy  per  unit  new  surface  generated  differed  by  a  factor 

of  2  within  a  given  hypothesis,  and  exhibited  an  average  value  of  2.25 

-3  -3  -3  2 

x  10  ,  about  6  x  10  and  approximately  11  x  10  j/mm  for  cases  (a), 

(b)  and  (c),  respectively.  The  comminution  efficiency ,  i.e.  the  frac¬ 
tion  of  initial  kinetic  energy  utilized  for  the  fragmentation  process, 
can  be  noted  in  Table  8  to  be  of  the  order  of  70  to  85  percent  and  to 

be  essentially  independent  of  initial  projectile  energy  for  the  range  of 
values  employed  in  the  present  investigation. 

f.  Relation  of  Results  of  Comminution  Theories 

A  characteristic  particle  size  for  the  ejecta  was  chosen  for 
each  test  as  the  mean  value  of  the  normal  distribution,  i.e.  y  ,  for 
purposes  of  comparison  with  the  comminution  theories  of  Kick,  Bond  and 
Rittinger.  The  characteristic  dimension  of  the  specimen,  to  which 
the  above  was  to  be  compared,  was  chosen  as  the  magnitude  of  the  vector 
resultant  of  the  crater  radius  and  depth.  Table  12  contains  the  particle 
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Hypothesis  (c),  all  cracks  are  pre-existing 


TABLE  12.  COMMINUTION  CONSTANTS  BASED  ON  QUASISTATIC  THEORIES 
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size,  energy  of  comminution  E 
’  OJ  com 

from  the  Kick,  Rittinger  and  Bond 

1/L2] 


C,  In 
k 


and  the  constants  C^,  0r  and  predicted 
equations,  respectively,  given  by 

(18)  Kick 


E  -  C  In  —  r^-1  (19)  Rittinger 

r  Ll2  liJ 

E  “  2Cb[v/fe  (1  }]  (20)  Bond 

where  and  are  the  characteristic  initial  and  final  dimensions 
(length,  area,  or  volume)  and  E  is  the  energy  of  comminution.  Although 
the  magnitudes  of  these  empirical  constants  have  been  determined  for  the 
present  series  of  experiments,  based  on  a  scale  involving  a  length 
dimension,  and  appear  to  increase  with  increasing  input  energy,  the 
intent  of  the  equations  (18)  to  (20)  was  to  utilize  ratios  of  either 
volume  or  surface  area  rather  than  length,  and  thus  a  difference  in  the 
numerical  results  involving  one  or  even  more  orders  of  magnitude  might 
occur.  Furthermore,  since  the  block  size  for  these  equations  was  presumed 
to  be  semi-infinite,  the  initial  value  of  the  constants  is  based  on  an 
artificially  selected  magnitude.  Since,  furthermore,  the  relations  were 
designed  for  applications  involving  crushing  and  grinding  of  a  large  mass 
into  many  smaller  parts  under  quasi-static  loading,  it  is  likely  that 
the  present  results  utilizing  these  empirical  relations  cannot  be  extra¬ 
polated  to  other  impact  situations. 

However,  it  was  found  that  the  mean  particle  size  tends  to 
decrease  with  increasing  impact  energy,  as  might  be  expected.  Although 
the  conservation  of  mass  dictates  that  this  size  dimension  must  approach 
a  lower  limit,  the  present  range  of  tests  is  insufficient  to  provide 
this  information. 

2 .  Shale 

a.  Ejecta 

Since  the  earlier  tests  at  comparable  impact  velocities  indi¬ 
cated  that  minimal  quantities  of  ejecta  were  produced  in  the  shale 
specimens,  only  a  few  of  the  experiments  conducted  during  the  second 
year  of  operation  weie  photographed,  with  a  typical  result  shewn  in 
Figure  18.  The  cloud  formation  is  distinctively  different,  forming  a 
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I’ig.  18  -  I’i  lm  History  of 
Typical  ShnLe  Test 

Run  S  40 


thin  sheaf  rather  than  a  relatively  dense  conical  shell.  It  was  found  in 
the  present  tests  that  both  the  velocity  of  the  fastest  particles  and 
the  total  mass  of  the  ejecta  were  of  the  order  of  1/10  that  of  the 
striking  sphere;  thus,  the  total  energy  of  the  ejecta  amounted  to  no 
more  than  approximately  0.1  percent  of  the  striking  energy.  In  all 
photographs,  the  rebound  trajectory  of  the  striker  was  obscured  by  the 
ejecta  cloud,  but  the  terminal  location  of  the  ball  in  the  catcher  tank 
clearly  indicated  a  negligible  rebound  velocity, 
b.  Dynamic  Penetration 

A  summary  of  the  penetration  results  due  to  the  impact  of  the 
steel  sphere  on  the  shale  specimens  is  presented  in  Table  13.  Recalling 
the  appearance  of  the  penetration  configuration  as  a  cylindrical  hole 
with  a  hemispherical  bottom,  with  either  very  small  or  no  dish  around 
the  entrance,  it  is  possible  to  attempt  to  model  the  compaction  pheno¬ 
menon  as  a  rigid  projectile  of  mass  in  entering  a  linearly  viscous 
medium  characterized  by  constant  c,  neglecting  elastic  compression 
rebound,  ejecta  and  wave  phenomena.  The  equation  for  this  simple  model 
is 


mx  +  cx  =  0  (21) 
where  x  is  the  penetration.  The  energy  absorbed  by  the  viscous  forces 
is  then  equal  to  the  initial  kinetic  energy  of  the  striker  and  is  given 


E  =  -jT  cx  dx  (22) 

where  P  is  the  maximum  penetration  depth.  The  equation  of  motion  has 
the  general  solution 


x  =  x  e_ct//ra  +  C  (23) 

o  o  v  ' 

where  the  constants  of  integration,  xq  and  Cq,  are  evaluated  from  the 

conditions 


t  =  0,  x  =  0,  x  =  vq  and  t  =  oo,x  =  P,x"0  (24) 

Upon  substitution  in  Equation  (23),  the  viscous  constant  c  can  be 
evaluated  as 


c  *  mv0/1? 


(25) 
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TABLE  13. 

DYNAMIC  SHALE  PENETRATION  DATA  FOR  A  *«-IN. 
SPHERE 

DIAMETER 

STEEL 

Shale (1)’ 

...  P,  Penetration 
Depth,  mm 

v.  Velocity 
m/s 

E,  Energy 
J 

,  P/mv 
mm/g-m/s 

P*/E 
mm  /J 

S20 

1.74 

51 

1.4 

.0341 

2.163 

S40 

3.53 

107 

6.0 

.0330 

2.077 

S60 

4.32 

131 

8.9 

.0330 

2.097 

S80 

5.08 

152 

12.0 

.0334 

2.151 

S100 

5.43 

169 

14.9 

.0316 

1.979 

S120 

5.68 

181 

16.9 

.0314 

1.909 

S140 

5.94 

194 

19.4 

.0306 

1.819 

S160 

6.53 

205 

22.2 

.0319 

1.921 

S180 

7.06 

213 

23.7 

.0331 

2.103 

S200 

7.23 

227 

26.7 

.0319 

1.958 

S220 

7.98 

235 

28.8 

.0340 

2.211 

S240 

8.12 

242 

30.4 

.0336 

2.161 

A9 

31.75 

852 

378 

.0373 

2.667 

AlO 

31.75 

881 

404 

.0360 

2.495 

All 

31.75 

570.59 

169.30 

.0535 

5.958 

A12 

88.90 

1718 

1535 

.0517 

5.149 

A13 

29.46 

512 

136 

.0575 

6.382 

A14 

12.45 

342 

153 

.0364 

1.013 

A15 

32.00 

896 

417 

.0357 

2.456 

A16 

29.21 

563 

165 

.0519 

5.171 

A17 

29.21 

567 

167 

.0515 

5.109 

A18 

26.42 

569 

168 

.0464 

4.155 

A19 

24.89 

554 

160 

.0449 

3.872 

A20 

29.97 

608.69 

192.66 

0.0473 

4.663 

A21 

25.40 

685 

244 

.0371 

2.644 

A22 

25.40 

591.01 

181.63 

0.0413 

3.552 

A23 

22.10 

497 

128 

.0445 

3.816 

A24 

38.35 

945.18 

464.55 

0.0390 

3.166 

75 


TABLE  13  (cont'd).  DYNAMIC  SHALE  PENETRATION  DATA  FOR  A  k-IN.  DIAMETER 

STEEL  SPHERE 


Shale 

P,  Penetration 
Depth,  mm 

v.  Velocity 
m/s 

E,  Energy, 

J 

P/mv 

mm/g-m/s 

P‘/E 
mm /J 

A25 

39.12 

966 

485 

.0405 

3.160 

A26 

62.23 

965.61 

484.85 

0.0620 

7.987 

A27 

54.36 

2487.17 

3216.72 

0.0210 

0.918 

A28 

60.45 

1501 

1172 

.0403 

3.118 

A29 

17.02 

537 

150 

.0317 

1.931 

A30 

16.00 

477 

118 

.0335 

2.169 

A31 

44.45 

939 

458 

.0473 

4.314 

A3  2 

47.42 

974 

493 

.0436 

4.561 

A3  3 

16.00 

527 

144 

.0304 

1.776 

A3  5 

14.99 

537 

150 

.0279 

1.498 

A36 

63.50 

1601 

1333 

.0397 

3.025 

A3  7 

63.50 

1618 

1361 

.0392 

2.963 

A38 

37.34 

843 

370 

.0443 

3.768 

A3  9 

38.35 

875 

398 

.0438 

3.695 

A40 

20.32 

502 

131 

.0405 

3.152 

A41 

20.07 

550 

157 

.0365 

2.566 

A42 

14.99 

449 

105 

.0334 

2.148 

A433 

64.008 

1603.55 

1337.12 

0.0384 

3.064 

A45 

78.99 

1609 

1346 

.0491 

4.612 

A46^ 

39.878 

874.78 

397.92 

0.0438 

3.996 

refers  to  tests  conducted  at  the  University  of  California,  Berkeley 

^A  refers  to  tests  conducted  at  the  Naval  Weapons  Center,  China  Lake, 
California. 
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Actual  integration  of  the  energy  equation  (22)  yields 

E  =  (1/2) mvQ2  (26) 

as  expected,  and  substitution  for  vq  yields 

(l/2)mvo2  =  Q  P2  where  Q  =  c2/2ra  i  (27) 

The  calculated  values  of  Q  and  c  for  each  test  are  reproduced  in  Table 

2 

13,  their  average  being  0.489  j/mm  and  30.3  kg/sec,  respectively,  with 
very  little  scatter.  It  is  expected  that  significant  deviations  will 
occur  in  this  model  when  the  initial  velocity  is  so  low  that  a  penetra¬ 
tion  depth  of  at  least  one  projectile  radius  is  not  obtained,  nor  under 
circumstances  of  high  initial  velocities  when  ejecta  formation,  gross 
cracking  of  the  sample,  and  even  backward  spall  are  in  evidence. 

In  Reference  [8],  a  post  mortem  section  of  a  5-1/2  inch  thick 
shale  sample  with  a  diameter  of  5-1/2  inch,  encased  in  plaster  of  Paris 
as  in  the  current  tests,  struck  by  a  1/4  inch  diameter  steel  sphere  at 
a  velocity  of  5300  ft/sec  was  exhibited;  the  hole  had  been  filled  with 
a  glue  prior  to  sectioning. In  this  photograph  reproduced  in  Figure  19, 
the  hole  appears  as  a  truncated  cone  with  a  small  cone  angle  and  a  dish 
at  the  entrance  region;  the  projectile  was  embedded  at  the  tip  of  this 
cone  whose  diameter  was  slightly  smaller  than  that  of  the  sphere,  indi¬ 
cating  some  recovery  of  the  compacted  material.  The  discolored  zone 
appearing  as  a  bell-like  shape  around  the  hole  is  not  necessarily  an 
indication  of  the  compacted  region,  as  it  may  represent  the  extent  of 
the  glue  penetration.  The  marked  difference  in  appearance  between  this 
photograph  and  the  shale  crater  produced  at  considerably  lower  velocities 
indicates  the  difference  in  the  penetration  patterns  occurring  in  differing 
velocity  regimes. 

c.  Quasi-static  Tests 

The  results  of  the  tests  on  the  Instron  machine  are  summarized 
in  Table  14.  Significant  scatter  and  no  consistency  with  respect  to  rate 
was  obtained  in  the  measured  energy  required  to  produce  a  given  inden¬ 
tation  under  identical  velocities  of  penetration;  however,  the  data 
indicates  a  drastically  greater  energy  requirement  for  a  given  depth  of 
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TABLE  14.  QUASISTATIC  PENETRATION  TESTS  ON  SHALE 
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TABLE  14  (con'd).  QUASISTATIC  PENETRATION  TESTS  ON  SHALE 
P,  Penetration  Depth,  mm 
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Fig.  19  Diametral  Section  of  a  5^  -in.  D.  x  -in.  Shale 
Core,  struck  by  a  ^  -in.  Diameter  Steel  Sphere  at 
a  Velocity  of  5300  ft/sec.  Crater  filled  with  a 
Catalyzed  Methyl  Methacrylate  Glue  Prior  to 
Sectioning. 
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indentation  in  corresponding  impact  tests  as  might  be  expected.  While 
in  the  initial  range  of  penetration  prior  to  a  depth  of  one  radius, 
there  appears  to  be  no  obvious  relationship  between  energy  and  penetra¬ 
tion  depth,  a  nearly  linear  equation  seems  to  exist  between  these  quan¬ 
tities  beyond  this  domain. 

d.  Correlation  of  Penetration  Data 

An  examination  of  the  penetration  phenomenon  in  shale  at  all 
impact  velocities  employed  in  the  present  study  ranging  from  51  to  2487 
m/sec  indicates  that  the  viscous  model  is  a  reasonable  approximation  of 
the  behavior  pattern,  whether  the  trajectory  was  normal  or  inclined  to 
the  target  surface.  The  data  from  the  rifled  and  unrifled  projectiles 
fired  from  a  tank  gun  are  not  included  here  due  to  sizable  inconsis¬ 
tencies  in  penetration  depths  relative  to  the  other  results  occasioned 
by  tumbling  that  could  not  occur  in  the  case  of  spherical  projectiles. 

The  curve  of  penetration  depth  vs.  initial  projectile  momentum  is  presented 
in  Figure  20  and  indicates  that  a  single  line  with  a  value  c  =  29.5  kg/ sec 
is  a  good  representation  for  all  but  the  highest  velocity  tests,  where 
three  anomalies  occur:  the  two  shots  with  penetration  depths  greater 
than  anticipated  involved  special  experimental  circumstances,  one  being 
a  splitting  of  the  sample  and  the  other  a  back-up  by  a  hard-steel  cylin¬ 
der.  The  highest  impact  velocity  employed  resulted  in  a  penetration 
depth  considerably  smaller  than  expected  on  the  basis  of  a  linear  repre¬ 
sentation  as  the  result  of  projectile  fragmentation.  The  initial  pro¬ 
jectile  energy  has  been  plotted  against  the  square  of  the  penetration 

depth,  as  shown  in  Figure  21;  this  diagram  indicates  that  while  a  single 

2 

line  with  constant  Q  =  0.339  J/mm  approximates  the  entire  span  of  the 

data,  a  better  description  of  the  experimental  results  is  given  by  one 

line  for  the  low-velocity  tests  and  another  for  the  higher  velocity 

2  2 

experiments,  with  values  of  Q  =  0.489  J/mm  and  Q  =  0.284  J/mm  ,  respec¬ 
tively.  The  overall  pattern  of  the  results  agrees  at  least  qualitatively 
with  that  of  Reference  [15]. 
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Fig.  21  -  Square  of  Penetration  Depth  as  a  Function  of 
the  Initial  Projectile  Energy  for  Shale 
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SECTION  V 
SUMMARY 

This  document  represents  the  final  technical  report  for  contract  No. 
DNA001-73-C-0226  between  the  University  of  Califqrnia,  Berkeley  and  the 
Defense  Nuclear  Agency  and,  together  with  report  DNA  3416Z,  describes  the 
accomplishments  achieved  during  the  two-year  program  of  the  project.  The 
work  consisted  of  laboratory  investigations  of  the  inpact  of  a  spherical  steel 
projectile  on  disks  of  diorite  and  shale  at  velocities  ranging  frcm  51  to  1718 
m/sec  and  a  selected  series  of  ordinance  firings  in  the  field.  The  field 
tests  and  higher  velocity  laboratory  tests  wer<~  performed  by  the  personnel  of 
the  Naval  Weapons  Center,  China  Lake,  and  were  described  completely  in  Refer¬ 
ence  [8].  On  the  basis  of  the  overall  results,  it  is  concluded  that: 

1.  The  comminution  theories  of  Bond,  Kick  and  Rittinger,  designed  for 
quasi-static  crushing  of  brittle  rock,  are  not  directly  applicable  to  projec¬ 
tile  impact  on  a  target  modeled  as  a  semi-infinite  solid.  Since  the  under¬ 
lying  assumption  of  these  laws  is  that  the  basic  shape  of  the  particle  remains 
unchanged,  it  is  obvious  that  any  process  which  leaves  all  but  a  relatively 
small  percentage  cf  the  parent  configuration  completely  intact  falls  outside 
the  scope  of  these  theories. 

2.  The  particle  size  distribution  obtained  from  the  present  comminution 
process  in  diorite  appears  to  follow  a  normal  logarithmic  pattern  permitting 
its  description  in  terms  of  the  mean  and  the  standard  deviation,  somewhat 
analogous  to  what  was  found  in  crushing  and  grinding  operations.  Higher 
impact  velocities  produced  smaller  mean  sizes  of  the  ejecta. 

3.  The  ballistic  pendulum  arrangement  originally  employed  could  be 
expediently  replaced  by  a  fixed  test  stand  since  the  energy  consumed  in  the 
swing  was  found  to  be  a  completely  negligible  fraction  of  the  initial  value. 

4.  Impact  on  different  types  of  rock  produced  different  shapes  of  the 
emitted  ejecta.  In  the  case  of  diorite,  this  formed  a  hollow  cone  whose 
kinetic  energy  was  found  to  be  about  10  percent  of  the  initial  value;  in  the 
case  of  limestone  and  shale,  a  rodlike  emission  was  observed,  but  with  a 
kinetic  energy  of  only  about  0.1  percent  of  the  impinging  striker. 
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5.  The  rebound  velocity  of  the  striker  amounted  to  about  2  percent  in 
the  case  of  diorite  and  a  negligible  fraction  in  the  case  of  the  softer 
materials.  Larger  ejected  fragments,  initially  separately  examined,  also 
failed  to  represent  more  than  a  vanishingly  small  portion  of  the  available 
energy. 

6.  The  energy  propagated  in  the  form  of  strain  waves  clustered  around 
a  value  of  about  10  percent  of  the  original  energy.  The  strain  wave  was  the 
agent  that  produced  additional  cracks  beyond  those  preexisting  outside  the 
domain  of  the  crater;  however,  the  area  associated  with  these  cracks  could 
only  be  estimated  in  terms  of  upper  and  lower  bounds. 

7.  The  energy  of  comminution  for  the  diorite  samples  was  found  to  be 
in  the  range  from  70  to  85  percent  of  the  total. 

8.  Although  enormous  scatter  was  observed  in  the  quasi-static  experi¬ 
ments  on  shale,  the  dynamic  tests  at  lower  velocities  conformed  very  well  to 
a  simple  linearly  viscous  model  that  was  also  a  reasonable  approximation  for 
the  higher  velocity  regimes. 

9.  No  permanent  deformation  of  the  projectiles  was  observed  in  any  of 
the  tests  of  the  present  program. 

10.  The  ejecta  surface  area  was  computed  on  the  basis  of  the  measured 
crater  mass  and  three  assumed  shapes  of  the  ejecta  characterized  by  a  single 
dimensional  parameter;  these  geometries  were  believed  to  be  a  reasonable 
approximation  of  the  captured  fragments  based  on  their  microscopic  analysis. 

The  original  proposal  addressed  itself  to  the  examination  of  comminution 
in  soft  rocks  by  various  experimental  techniques  that  would  provide  separately 
information  concerning  the  diverse  phenomena  involved  in  the  overall  process. 
Both  shale  and  limestone  samples  were  examined;  however,  when  it  was  noted 
that  comminution  did  nor  play  a  major  role  in  the  laboratory  impact  tests  on 
soft  rock,  an  expanded  sequence  of  experiments  was  undertaken  for  the  hard 
rock  diorite  at  these  speeds  in  order  to  study  cases  where  ejecta  formation 
appeared  to  be  important.  Nevertheless,  in  tests  on  shale  involving  a  range 
of  initial  kinetic  energies  from  1.35  J  to  3.65  x  10  6  J,  the  penetration 
depth  scaled  reasonably  well  with  both  initial  momentum  and  initial  energy. 
Particle  size  distributions  have  been  obtained  both  in  the  lower  and  the 
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higher  velocity  regions  and  recapitulated  in  Reference  [8].  Thus,  results 
have  been  obtained  in  all  areas  stated  in  the  proposal,  and  estimates  of  the 
energy  required  to  produce  unit  new  surface  have  been  advanced. 

It  was  originally  indicated  that  the  period  of  investigation  was  too 
short  to  permit  a  caiprehensiw  analysis  and  modeling  of  the  inpact  and  com¬ 
minution  phenomenon  and  that  the  present  investigation  could  serve  at  most 
as  a  point  of  departure  far  such  an  all-inclusive  program,  and  this  turned 
out  to  be  the  case.  Severe  experimental  difficulties  were  encountered  in  the 
collection  and  sizing  of  all  ejecta,  whether  at  the  lower  or  higher  inpact 
velocities.  The  determination  of  the  new  surfaces  formed  by  the  cracks 
generated  by  the  passage  of  the  transient  produced  by  the  impact  could  not  be 
accomplished  with  any  degree  of  accuracy  or  consistency;  this  is  believed  to 
be  due  to  the  inhomogeneous  distribution  of  the  preexisting  microcracks. 
Finally,  further  information  is  needed  for  the  construction  of  a  proper  me¬ 
chanical  model  that  will  predict  the  microscopic  phenomena  involved  in  the 
cracking  fracture  and  ejecta  formation  of  hard  rocks  such  as  diorite.  In 
the  case  of  shale,  where  oonpaction  appears  to  be  an  additional  aspect  of 
the  contact  problem,  it  is  possible  that,  at  least  in  the  lower  ranges  of 
inpact  velocity,  a  simple  model  .based  on  viscous  considerations  could  be  re¬ 
fined  to  provide  adequate  predictive  capabilities. 
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Appendix  A  -  Computer  Program  for  Crack  Surface  Area 
Evaluation 
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Sandla  Laboratories 

ATTN:  Doc.  Con.  for  3141,  Sandia  Rpt.  Coll. 

U.  S.  Energy  Research  t  Development  Administration 

Albuquerque  Operations  Office 

ATTN:  Doc.  Con.  for  Tech.  Lib. 

U.S.  Energy  Research  t  Development  Administration 

Division  of  Headquarters  Services 

ATTN:  Doc.  Con.  for  Class.  Tech.  Lib. 

U.S.  Energy  Research  b  Development  Administration 

Nevada  Operations  Office 

ATTN:  Doc.  Con.  for  Tech.  Lib. 

Union  Carbide  Corporation 

ATTN:  Civil  Def.  Res.  ProJ. 

ATTN:  Doc.  Con.  for  Tech.  Lib. 
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OTHER  GOVERNMENT  AGENCIES 


DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 


Bureau  of  Mines 

ATTN:  Tech.  Lib. 

Department  of  the  Interior 

Bureau  of  Mines 

ATTN:  Tech.  Lib. 

Department  of  the  Interior 

U.S.  Geological  9urvey 
ATTN:  J.  H.  Healy 
ATTN:  Cecil  B.  Raleigh 

DEPARTMENT  OF  DEFENSE  CONTRACTORS 

Aerospace  Corporation 

ATTN:  Tech.  Info.  Services 

Agbabian  Associates 

ATTN:  M.  Agbabian 

Applied  Theory,  Inc. 

2  cy  ATTN:  John  G.  Trullo 

Avco  Research  &  Systems  Group 

ATTN:  Research  Library,  A-830,  Rm.  7201 

Battelle  Memorial  Institute 
ATTN:  Tech.  Lib. 

The  BDM  Corporation 
ATTN:  Tech.  Lib. 

The  Boeing  Company 

ATTN:  Aerospace  Library 

California  Research  6  Technology,  Inc. 

ATTN:  Ken  Kreyenhagen 
ATTN:  Tech.  Lib. 

University  of  California 
ATTN:  G.  Sacl.man 
50  cy  ATTN:  W.  Goldsmith 

Calspan  Corporation 
ATTN:  Tech.  Lib. 

Civil /Nuclear  Systems  Corporation 
ATTN:  Robert  Crawford 

University  of  Dayton 

ATTN:  Hal  lock  F.  Swift 

University  of  Denver 

Colorado  Seminary 

ATTN:  Sec.  Officer  for  Tech.  Lib. 

ATTN:  Sec.  Officer  for  J.  Wlsotskl 

EG&G,  Inc. 

Albuquerque  Division 
ATTN:  Tech.  Lib. 

General  American  Transportation  Corporation 
ATTN:  G.  L.  Neidhardt 

General  Electric  Company 

TEMPO-Center  for  Advanced  Studies 
ATTN:  DASIAC 


IIT  Research  Institute 

ATTN:  Milton  R.  Johnson 
ATTN:  R.  E.  Welch 
ATTN:  Tech.  Lib. 

Institute  for  Defense  Analyses 

ATTN:  IDA,  Librarian,  Ruth  S.  Smith 

Kaman  AviDyne 

Division  of  Kaman  Sciences  Corporation 
ATTN:  E.  S.  Crlscione 
ATTN:  Norman  P.  Hobbs 
ATTN:  Tech.  Lib. 

Kaman  Sciences  Corporation 
ATTN:  Library 

Lockheed  Missiles  &  Space  Company 

ATTN:  Tom  Geers,  Dept.  52-33,  Bldg.  205 
ATTN:  Tech.  Info.Ctr.,  D/Coll. 

Lovelace  Foundation  for  Medical  Education 

ATTN:  Asst.  Dir.  of  Res.,  Robert  K.  Jones 
ATTN:  Tech.  Lib. 

McDonnell  Douglas  Corporation 
ATTN:  Robert  W.  Halprin 

Merritt  Cases,  Incorporated 
ATTN:  J.  L.  Merritt 
ATTN:  Tech.  Lib. 

The  Mitre  Corporation 
ATTN:  Library 

Nathan  M.  Newmark 

Consulting  Engineering  Services 
ATTN:  Nathan  M.  Newmark 

Physics  International  Company 

ATTN:  Doc.  Con.  for  Larry  A.  Behrmann 
ATTN:  Doc.  Con.  for  Tech.  Lib. 

ATTN:  Doc.  Con.  for  Dennis  Orphal 
ATTN:  Doc.  Con.  for  Robert  Swift 
ATTN:  Doc.  Con.  forE.  T.  Moore 
ATTN:  Doc.  Con.  for  Charles  Godfrey 
ATTN:  Doc.  Con.  for  Fred  M.  Sauer 

R  &  D  Associates 

ATTN:  William  B.  Wright,  Jr. 

ATTN:  Jerry  Carpenter 
ATTN:  J.  G.  Lewis 
ATTN.  Tech.  Lib. 

ATTN:  Albert  L.  Latter 
ATTN:  Sheldon  Schuster 
ATTN:  Henry  Cooper 
ATTN:  Harold  L.  Brode 
ATTN:  Bruce  Hartenbaum 

Science  Applications,  Inc. 

ATTN:  D.E.  Maxwell 
ATTN:  David  Bernstein 

Science  Applications,  Inc. 

ATTN:  R.  Seebaugh 
ATTN:  John  Mansfield 
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DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 

Science  Applications,  Inc. 

ATTN:  Tech.  Lib. 

ATTN:  Michael  McKay 

Science  Applications,  Inc. 

ATTN:  R.  A.  Shunk 

Shock  Hydrodynamics,  Inc. 

A  Div.  of  Whittaker  Corporation 
ATTN:  L.  Zernow 

Southwest  Research  Institute 
ATTN:  A.  B.  Wenzel 
ATTN:  Wilfred  E.  Baker 

Stanford  Research  Institute 
ATTN:  Carl  Petersen 
ATTN:  George  R.  Abrahamson 
ATTN:  Burt  R.  Gasten 

Systems,  Science  &  Software,  Inc. 

ATTN:  Robert  T.  Allen 
ATTN:  Donald  R.  Grine 
ATTN:  Tech.  Lib. 

ATTN:  Ted  Cherry 
ATTN:  Thomas  D.  Rlney 

Terra  Tek,  Inc. 

ATTN:  Sidney  Green 
ATTN:  Tech.  Lib. 

Tetra  Tech,  Inc . 

ATTN:  LI- San  Hwang 
ATTN:  Tech.  Lib. 

TRW  Systems  Group 

ATTN:  Benjamin  Sussholtz 
ATTN:  Paul  Lieberman 
ATTN:  William  Rowan 
ATTN:  Norm  LIpner 
ATTN:  Tech.  Info.  Ctr.,  S-1930 
ATTN:  Pravin  Bhntta 
ATTN:  Jack  Farrell 


DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 

TRW  Systems  Group 
San  Bernardino  Operations 
ATTN:  Greg  Hulcher 

Universal  Analytics,  Inc. 

ATTN:  E.  I.  Field 

The  Eric  H.  War.g 

Civil  Engineering  Research  Facility 
ATTN:  Neal  Baum 
ATTN:  Larry  Bickle 

Washington  State  University 
Administrative  Office 

ATTN:  Arthur  Miles  Hohorf  for  George  Duval 

Weldltnger  Assoc.  Consulting  Engineers 
ATTN:  Melvin  L.  Baron 

Weldltnger  Assoc.  Consulting  Engineers 
ATTN:  J.  Isenberg 

Westinghouse  Electric  Company 
Marine  Division 

ATTN:  W.  A.  Volz 
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